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Abstract 

We investigate well-posedness for martingale solutions of stochastic differential equa¬ 
tions, under low regularity assumptions on their coefficients, widely extending the results 
first obtained by A. Figalli in [Fig08]. Our main results are a very general equivalence 
between different descriptions for multidimensional diffusion processes, such as Fokker- 
Planck equations and martingale problems, under minimal regularity and integrability 
assumptions, and new existence and uniqueness results for diffusions having weakly dif¬ 
ferentiable coefficients, by means of energy estimates and commutator inequalities. Our 
approach relies upon techniques recently developed jointly with L. Ambrosio in [AT14], 
to address well-posedness for ordinary differential equations in metric measure spaces: 
in particular, we employ in a systematic way new representations and inequalities for 
commutators between smoothing operators and diffusion generators. 


1 Introduction 

Aim of this article is to study well-posedness (i.e., existence, uniqueness and stability) for 
martingale solutions of stochastic differential equations 

dXt = b{t,Xt)dt + a{t,Xt)dWt, t€(0,r), (1) 

providing in particular new results, under low regularity assumptions on the coefficients b : 
(0,T) a: (0,T) x M"* ^ 

The classical subject of martingale problems dates back at least to [SV06], where it was 
hrst shown that continuous and uniformly elliptic covariances a = acr*’s allow for uniqueness 
results which have no counterpart in the usual (It6-)Cauchy-Lipschitz theory, provided that 
the solution to (1) is understood in a sufficiently weak sense. Since then, the theory has been 
growing, due to its robustness and strong connections with the theory of semigroups and 
parabolic PDE’s, also in abstract (metric) frameworks, see e.g. [EK86]. 

Our primary goal here is to show that the techniques originally developed in [AT14] can 
be extended to the stochastic theory as well as specialized to the Euclidean setting, to extend 
in a systematic way the results established in the seminal paper [Eig08]. Actually, most of 
such techniques, tailored to study well-posedness problems for ordinary differential equations 
in metric measure spaces (possibly infinite-dimensional) are also well-suited also to the study 
of diffusions in metric measure spaces, as developed in the author’s PhD dissertation [Trel4]. 
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However, in this paper, we deal uniquely with Euclidean spaces: among various motivations, 
besides that a wider audience could be mainly interested in this setting, this allows us to 
compare new results and techniques with alternative approaches. Finally, Euclidean spaces 
are a useful “intermediate” step for the infinite dimensional theory, e.g. by cylindrical ap¬ 
proximations; the theory developed here is also instrumental to the developments in [TrelS]. 

Therefore, in this article, we adopt the same point of view as in [Fig08], where precise 
connections between well-posedness of PDF’s and martingale problems are settled, in par¬ 
ticular for a wide class of diffusion having not necessarily continuous nor elliptic coefficients, 
provided that some Sobolev regularity holds. Of course, well-posedness has to be understood 
“in average” with respect to ..Sf'^-a.e. initial condition (here and below, is Lebesgue mea¬ 
sure on M'’*). More precisely, a formalization akin to that of DiPerna-Lions (see e.g. [AC14] 
for an account of the deterministic theory) is introduced, the main objects being Stochastic 
Lagrangian Flows, i.e., Borel families {'n{x))xeRd of probability measures on (^([O,T];M^), 
such that 

(i) r]{x) solves (1), starting from x ai t = 0, for x G 

(ii) the push-forward measures (et){j f r/(x) dJf^^lx), where et is the evaluation map at t G 
[0,T], are absolutely continuous with respect to with uniformly bounded densities. 

Let us stress the fact that, as in the deterministic theory, uniqueness is understood for flows, 
thus in a selection sense: we are not claiming well-posedness for ..Sf'^-a.e. initial datum. More¬ 
over, we remark that, although the conditions above might read as perfect analogues of the 
notion of Regular Lagrangian flows [AC14, Definition 13], Stochastic Lagrangian Flows are 
not necessarily (neither expected to be) deterministic maps of the initial point only; this is 
evident when a = 0 above and any probability concentrated on possibly non-unique solu¬ 
tions to the ODE give rise to a solution to the martingale problem. Despite this discrepancy, 
such a theory provides rather efficient tools to study stochastic differential equations under 
low regularity assumptions, in Euclidean spaces, and, together with [LBL08], which deals 
with analogous issues from a PDF point of view, has become the starting point for further 
developments, among which we quote [RZ10,Luol3,FLT10,Zhal3]. 

Before we proceed with a more detailed description of our results and techniques, let us 
stress the fact that we are concerned uniquely with martingale problems, so we do not address 
nor compare our results with those obtained for strong solutions of equations under low regu¬ 
larity assumptions on the coefficients (see the seminal paper [VerSO] and [KR05,DFPR13] for 
more recent results). Rigorous correspondences between martingale (or weak) and strong so¬ 
lutions may be provided by the classical Yamada-Watanabe theorem [YW71] (and extensions, 
see e.g. [Kur07]). Moreover, the literature on Fokker-Planck equations for general measures 
is so vast that we must limit ourselves to a comparison of our results only with those which 
are strongly related and look similar in techniques and mathematical contents: this is done 
in Section 3.3. 

We proceed with a brief description of our contributions developed below, which can be 
split into two parts, roughly corresponding to Section 2 (toghether with Appendix A) and 
Section 3. 

In the first part, we investigate the problem of abstract equivalence between “Eulerian” 
and “Lagrangian” descriptions for multidimensional diffusion processes, where by the former 
we mean by Fokker-Planck equations and the latter consists of solutions to martingale prob¬ 
lems. Although such a correspondence can not be considered novel and many ideas can be 
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traced back at least to [Amb04] in the theory of ODE’s and DiPerna-Lions flows, as well 
as [KS98] for cadlag martingale problems, to our knowledge, here we provide for the first time 
general results, under somewhat minimal integrability assumptions on coefficients as well as 
on solutions. Moreover, we choose to state and prove our results in such a way that they 
can be translated with a minimal effort to the case of general metric measure spaces, that we 
address in [Trel5]. 

In this part, the crucial result is Theorem 2.5, which provides a so-called “superposition 
principle”, i.e., a (non-canonical) way to lift any probability-valued solution of a Fokker-Planck 
equation to some solution of the corresponding martingale problem. Here, “to lift” means that 
the 1-marginals of the process which solve the martingale problem coincide with the given 
solution of the Fokker-Planck equation. Results in a similar spirit appear quite often in the 
literature (see also the comments just below the statement of Theorem 2.5) and could be traced 
back to L.C. Young’s theory of generalized curves. Technically, one could start from already 
known results such as [Fig08, Theorem 2.6] or [KS98, Theorem 4.9.17] to provide a slightly 
shorter proof, but we preferred to postpone an almost self-contained derivation in Appendix 
A: indeed, even if we rely on the results quoted above, it turns out that one has to settle 
non-trivial technical problems. In particular, an underlying result is Theorem A.2, where 
we establish an estimate for the modulus of continuity of solutions to martingale problems 
under somewhat minimal integrability assumptions (based on a refined Levy-type estimate); 
an alternative but less effective approach, based on fractional Sobolev spaces, was developed 
in [Trel4]. Finally, we point out that we exploit a technique originally developed in [AT14, 
Theorem 7.1], in case of cylindrical approximations, to move from bounded coefficients to 
possibly unbounded ones. 

In the second part, we address the problem of well-posedness for Fokker-Planck equations, 
providing sufficient conditions assuming Sobolev regularity of the coefficients. We mainly fo¬ 
cus on uniqueness issues, which are settled by means of energy or estimates, formally 
satisfied by any weak solution, under suitable bounds on the divergence of the driving coeffi¬ 
cients; such an approach could be hardly considered novel, as it was already present in [DL89], 
for transport equations. However, our main contribution consists in a novel and systematic 
approach to the estimate of the error terms arising in the approximation procedure, to obtain 
so-called commutator inequalities; see Section 3.4 for a brief account of the method as well as 
complete proofs of our crucial resuts. It turns out that, essentially by means of the same tech¬ 
nique, we are able to deal with Sobolev derivations (Lemma 3.4), Sobolev diffusions (Lemma 
3.5) as well as with time-dependent elliptic diffusions (Lemma 3.6). Such a technique, which 
ultimately consists in choosing a Markov semigroup as a smoothing operator and relying on 
duality arguments as well as an interpolation a la Bakry-Emery, has also the advantage of 
being completely “Eulerian” and “coordinate free”. Let us point out also in this case that 
it was first developed in [AT14] to deal with an analogue problem for derivations in metric 
measure spaces. 

In conclusion, we state and prove two well-posedness results; Theorem 3.1, for diffusions 
(1) having possibly degenerate coefficients, assuming first order Sobolev regularity for the drift 
b and second order Sobolev regularity for the infinitesimal covariance a = aa* (together with 
uniform bounds on their divergence); Theorem 3.1, for the bounded elliptic case, i.e. A < 
a{v,v) < A \v^ for every v G with t ^ at Lipschitz, where (roughly speaking) regularity 
assumptions can be reduced of one order (i.e., no assumption on b, and hrst order Sobolev 
regularity for a). We regard such results as chief examples of the strength and versatility 
of our techniques for commutator estimates, and we point out that other interesting results 
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could arise in different situations, such as perturbations of elliptic generators which enjoy 
some ultra- (or hyper-) contractivity features, as well as the case of BV-regular coefficients, 
that we do not address here. 
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2 Diffusion processes and their equivalent descriptions 

In this section, we study abstract correspondences between “Eulerian” and “Lagrangian” de¬ 
scriptions for multidimensional diffusion processes, in particular with respect to well-posedness 
results. The main ideas involved are not entirely novel, but they widely extend those 
from [Fig08]: here we obtain results under minimal regularity and integrability assumptions. 
As already remarked in the introduction, on a technical side, a crucial tool is the superposition 
principle for diffusions. Theorem 2.5, whose proof is deferred to Appendix A. 

In Section 2.1, we introduce diffusion operators in W^, Fokker-Planck equations, martingale 
problems and flows; in Section 2.2 we study their equivalences. 


2.1 Definitions and basic facts 

Throughout, we use the following notation, for v, w ^ {d > 1) and A, B ^ 
d 

V ■ w = |u|^ = u-u, {v ® wY’^ ■.= v'^w\ for i, j G {1,... d}, 

i=l 


A: B = A^’^ B^’Y 1^1^ = A : A, A{v, w) = A : (v 0 w), 

i,j=^ 

and the following notation for differential calculus on (0,T) x R'^ (T G [0,oo)): 

/*(•) =/(b-), = = fortG[0,r],i,jG{l,., 


, 4 , 


d d 

V/ = VV = thus b • V/ = J2mf^nd a:V^f= 

i=l 

as well as the notation for Lebesgue measure on and V* for the distributional adjoint 
of V (i.e., V*6 = — div6 on vector fields). 

We write ^ (R*^) for the space of signed (real-valued) Borel measures on R*^ (with finite 
total variation), ./#■*■ (R'^) C ^ (R*^) for the cone of finite non-negative measures and C 

.y#^(R'^) for the convex set of Borel probability measures on R*^. We say that a curve u = 
i^t)t£{o,T) ^ .^(R'^) is Borel if, for every Borel set A C R*^, the curve t e-)■ t't(A) is Borel; we let 
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1^1 = (l^i|)iG(o,r) be the curve of total variation measures. A curve v = {i't)t&(o,T) ^ 
is narrowly continuous if, for every / G C'b(M'^), t >-)■ J fdr]t is continuous. 

Most of the quantities that we consider below are integrated with respect to the variable 
t G (0,T), with respect to when v = {i't)te{o,T) ^ is a Borel curve, we write 

\v\ dt for the Borel measure on (0,T) x A i-G \iy\ (A) = {At)dt, for A C (0,T) x R*^ 

Borel. For p, q G [l,oo] and a Borel curve = {Wt\)te{o,T) ^ ./^'''(R'^), the space is 

naturally defined and endowed with the Banach norm 


\\f\\LPL%iu) 


ll/(i,a:)|| 


L%{ut) 


< oo, 

Lf {dt) 


On the space (^([O, T]; R*^) (naturally endowed with the sup norm and its Borel cr-algebra), 
we let et : 7 i-T- := 7 (t) G R'^ be the evaluation map at f G [0,T]. The natural filtration on 
(^([O, T]; R'^) is the increasing family of cr-algebras T = (Tt)tg[o,T]) with Tt := a{es : s G [0,t]). 
Given ij G .^((^([O,T];R'^)), we always let pt ■= be the 1-marginal law at t G [0,T]. 

Notice that the family p := {pt)t£[o,T] ^ <^(R'^) is narrowly continuous. 

We let throughout £/ = C'^’^((0,T) x R*^)) (respectively, = C'b^((0,T) x R'^))) be 
the space of uniformly bounded (respectively, compactly supported) and continuously differ¬ 
entiable functions, once with respect to t G (0,T) and twice with respect to x G R'^, with 
uniformly bounded derivatives (as usual, the superscript (1,2) counts the number of deriva¬ 
tives with respect to {t,x), other superscripts may appear, with natural meaning). We prefer 
the “abstract” notation £/ and large parts of the theory can be developed when “test” func¬ 
tions are replaced by other classes (e.g. as developed throughout the monograph [KS98]). We 
endow £/ with the norm 


11/11^1,2= sup {\f{t,x)\ + \dtf{t,x)\ + \Vf{t,x)\ + \V‘^f{t,x)\} . 
7,a;)e(0,T)xR‘i 

Notice that, by uniform continuity, any f G extends to [0, T] x R*^. 

Throughout, we always let 

a = (a*b)f^^.^^ : (0,r) x R"* ^ Sym+(R'^), b = (6*)ti : (0,T) x R"* ^ R"*, (2) 

be Borel, where Sym_|_(R'^) is the space of symmetric, non-negative definite d x d matrices. 

We define diffusion operators in R*^, measure-valued weak solutions to Fokker-Planck 
equations and martingale problems on the interval [0,T]. Most of these notions are classical 
(for a brief historical account, see e.g. the introduction of [SV06]): for the sake of clarity we 
provide the definitions and prove some simple facts. 

Definition 2.1 (diffusion operator). We let L{= L{a,b)) be the linear differential operator 

^ T/ = ia:V2/ + 6-V/, 

with values in Borel maps on (0,T) x R^. 

We write Ltf := {Lf)t, for t G (0, T). As usual, the coefficients 6, a are referred as the 
drift of L and the infinitesimal covariance of T. If a = 0, then L reduces to a linear first-order 
operator, i.e. a derivation, and we say that we are in the deterministic case. 
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Given a diffusion operator L, we let the “Eulerian” description of evolution of particles 
“driven” by L consist of weak solutions of Fokker-Planck (or forward Kolmogorov) equations, 
in duality with Although our main interest lies in solutions to FPE’s that are narrowly 
continuous curves of probability measures, we introduce more general measure valued solu¬ 
tions, as they are useful, e.g. the space of solutions becomes linear. 

Definition 2.2 (weak solutions of FPE’s). A Borel curve v = {r't)te{o,T) ^ is a weak 

solution of the Fokker-Planck equation (FPE) 

dtvt = L^vt, on (0,r) X M"*, (3) 

if it holds 

J J {\at\ + \bt\) d\iyt\dt < oo (4) 

and, for every f G it holds 



[dtf{t,x) Lf{t,x)] dut{x) dt = 0. 


(5) 


With the notation introduced above, condition (4) can be restated as a, 6 G Ll,^{\v\). In 
what follows, we frequently omit to specify the operator £, that we regard as fixed. 


Remark 2.3. A density argument akin to [AGS08, Lemma 8.1.2] allows for proving that any 
solution V = (r't)tg(o,T) ^ ^(M'^) to (3) admits a unique narrowly continuous representative 
h = (P)tG[o,T]; with ut = ht, for jSf^-a.e. t G (0,T). Thanks to this fact, we may also say that 
the solution v starts from vq (or that vq is the initial law of v). Moreover, for every f G £/, 
it holds 


/ 


ftidht^ 



[dtf -h Ltf] dvtdt, for ti, t 2 G [0,T], ti < t 2 . 


( 6 ) 


Actually, one first proves that (6) holds for / G and then extends by density to / G . 
Since this last step requires the introduction of useful cut-off functions, we sketch it here, for 
later use. For ii > 1, we fix XR '■ [0; 1]; a smooth function with xr{x) = 1, for |x| < ii, 

Xr{x) = 0, for |x| > 2R, such that |Vy;_R| < and |V^Xi?| < 4i?“^. Given / G . 2 /, we let 

fR = fXR £ for which we assume that (6) holds. The chain rule entails 


^tfR = {Ltf )XR + ft^tXR + at{Vft,VxR), for t G (0,T), 


hence the bound 

\LtfR\ < \Ltf\ + \ft\\LtXR\ + \at\\Vft\ |Vxij| < C \\ft\U [|ai| + M . 

0 

Letting R —)■ 00 , by dominated convergence, we extend the validity of (6). 

Next, we introduce solutions of the martingale problem, following [SV06, Chapter 6]. In 
particular, we argue directly on the “canonical” space n = (^([O, T]; M'^), endowed with the 
evaluation process et(x) = 7{t), t £ [0)^”], and its natural filtration. 
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Definition 2.4 (solution of HP’s). A probability rj G is a solution of the 

martingale problem (MP) (assoeiated to L) if it holds 


{\bt\ oet + \at\ o et) dt 


1.70 


dr] < oo 


and, for every f G , the proeess 

[0,T] 3 t fto et - [ [dtfs + 'Csf] o Gsds 

Jo 


(7) 

( 8 ) 


is a martingale with respect to the natural filtration on (^([O, T]; R'^). 

Recall the notation pt = (et)#'?? £ t G [0,T], thus rjo is the initial law of rj. As for 

FPE’s, we usually omit to specify L, regarded as fixed. Let us remark that a density argument 
shows that it makes no difference to require that ( 8 ) is a martingale only for / G £/c- 

For any solution rj of the MP, the integrability assumption (7), which is equivalent to a, 
b G L^ig), entails that the process [0, T] 3 t Jp [dtfs + Ji^sf] ° ^sds is well defined, up to a 
T 7 -neghgible set, as continuous and progressively measurable process. In particular, it belongs 
to L^Jj], (3't)t), i.e. there exists an increasing sequence of stopping times r^, r]-a..s. converging 
towards T, such that [dtfs + L^/] o Cgds G L°°{r]), for every n > 1: it is sufficient to let 

Tn := F A inf G [ 0 , T] : J o + |as| o e^) ds > n| . 


We prefer throughout not to enlarge the filtration S' with r] negligible sets. This causes 
virtually no harm in the exposition, e.g. a martingale M = {Mt)t£[o,T] must be understood in 
the sense that it holds E[Mt|3's] = Mg, ij-a.s. for every s < t (and the 77 -negligible set could 
not belong to T^)- 

When a = 0, solutions to the MP reduce to probability measures concentrated on abso¬ 
lutely continuous solutions to the ordinary differential equation 

4:7t = hilt), for t G (0,r). 

dt 

Indeed, arguing as in [Fig08, Lemma 3.8], it turns out that the martingale ( 8 ) is constant. 
More generally, the quadratic variation process of ( 8 ) is t as{Vfg,'^fs)ds: this plays 

a crucial role in estimates for the modulus of continuity of the canonical process, see e.g. 
Corollary A.4. 

By integration of ( 8 ) with respect to 77 (i.e., taking expectation) we deduce that any 
solution 77 of the MP induces, by means of its 1-marginals ir]t)t&(o,T) ^ narrowly continuous 
solution of the FPE (3). A converse statement is provided by the following theorem, whose 
proof is deferred in Appendix A; in the next section, it plays a crucial role to connect various 
well-posedness results for FPE’s and MP’s. 

Theorem 2.5 (superposition principle). Let v = {vt)t£:{o,T] ^ ^(R'^) be a narrowly eontin- 
uous solution of (3). Then, there exists 77 whieh is a solution to the MP (assoeiated to the 
same diffusion operator L) such that, for every t G [0,r], it holds r]t = 
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In what follows, we refer to r] above as a superposition solution for u. 

We refer to this result as the superposition principle for diffusions; the terminology origi¬ 
nates in the deterministic literature of ODE’s, see [Amb04]: the solution rj can be non-trivially 
distributed among the possibly non-unique solutions to the ODE, thus introducing some “ran¬ 
domness” in an otherwise deterministic setting; these probability measures are nevertheless 
superpositions of deterministic paths. In the setting of diffusion operators, solutions are al¬ 
ready expected to be random, thus the term is justified only by extension, although it would 
be interesting, at least in some cases, to be able to distinguish between the two “sources of 
randomness”: this would require us to introduce concepts such as strong and weak solutions. 

As remarked in the introduction, Theorem 2.5 is a quite general result, only the integra- 
bility condition (4) being required, which is some sense minimal to give sense to FPE’s and 
MP’s (although one may slightly relax it by dealing with local martingale problems). Our re¬ 
sult extends [Fig08, Theorem 2.6], where only uniformly bounded coefficients are considered; 
let us mention that results in a similar spirit - that of L.C. Young’s theory of generalized 
curves - appear quite often in the literature, e.g. Echeverria’s theorem [EK86, Theorem 4.9.17] 
(see [KS98] for extensions) in the framework of martingale problems in spaces of cadlag paths, 
or Smirnov’s decomposition of 1-currents [Smi93] (see also [PS12] for an alternative approach, 
valid also in the case of metric currents). Our strategy of proof extends that of [Fig08, Theo¬ 
rem 2.6] and should be regarded as a (non-trivial) counterpart of [AGS08, §8.1 and §8.2] in the 
setting of multi-dimensional diffusions: although rather natural, the derivation is not imme¬ 
diate from the available literature (both from deterministic and stochastic), due to non-trivial 
technical points. The major difficulty in our proof is to provide estimates for the modulus of 
continuity of the canonical process (a problem that would appear also if we wanted to deduce 
it from Echeverria’s theorem). 

Next, we investigate some stability properties enjoyed by solutions of MP’s and FPE’s, 
with respect to suitable operations: their proofs are straightforward, so we omit them. 

Clearly, all the definitions above can be given with respect to any interval [S,T] in place 
of [0,T] (when it is not mentioned, we always refer to the interval [0,T]): solutions are then 
well-behaved with respect to the natural restriction map 

C([0,T];M'') 9 7^ = (7t)ie[5,T] G C([5,T];R'^). 

Proposition 2.6. Let 5, T G R, with 0 < S < T, and let rj G .^((^([O, T]; R'^)) he a solution 
of the MP. Let p : C'([5,T];R'^) —)■ [0,oo) be a uniformly bounded probability density (with 
respect to r]), measurable with respect to S's- 

Then, the push-forward {p'n)\[s,T] •= {\[s,T]) G ,^(C'([5,T];R'^) is a solution to the 
MP associated to L on[S,T]. 

The analogous property for FPE’s is obvious: if is a solution of (3), its restriction 

{r't)t£{s,T) is a solution of the FPE on (5, T) x R'’*. 

Solutions of FPE’s and MP’s are clearly stable with respect to convex combinations, as a 
consequence of Fubini’s theorem. 

Proposition 2.7. Let {Z,A,P) be a probability space and let {ri^)z£z C '^(C'[0,T];R‘’') be a 
Borel family, such that, for u-a.e. z G Z, ij,, is a solution of the MP (associated to a fixed 
diffusion operator L). Moreover, let 



(16tl -I- JatJ) dr}zdtdu{z) < oo 


(9) 



hold. Then, A i—)■ r]{A) = f r],,{A) du{z) is a solution of the MP (associated to L). 

A somewhat converse result, for disintegration with respect to the initial law, is a conse¬ 
quence of stability of martingales under conditional expectations with respect to the cr-algebra 
3 - 0 . 

Proposition 2.8. Let ij be a solution of the MP and let {il{x))x£Rd be a regular eonditional 
probability for rj with respeet to cq. Then, for Tjo-a.e. x G ri{x) is a solution of the MP 
assoeiated to L, with initial law 6x. 

We conclude this section by introducing a suitable notion of flow associated to a diffusion 
operator, roughly consisting of Borel families of solutions of the MP, for a (large, in some 
sense) set of initial conditions in Our aim is to study flows in the DiPerna-Lions sense 
(as extended by Figalli to MP’s), thus, we introduce the concept of “regular flow”, where 
regularity is usually some growth and/or absolute continuity condition on the 1-marginals, 
providing a selection criterion, yielding uniqueness in otherwise ill-posed problems. To study 
this notion in sufficient generality, we formulate such regularity conditions in terms of some 
set := f^[o,T] of narrowly continuous (probability curves that are) solutions of (3), which 
describe the “admissible” class of dynamics. With this notation, we refer to any n £ !R as a 
IR-regular solution of (3), and we say that solution to the MP is IR-regular if the curve of its 
1-marginals is a IR-regular solution of (3). We also let ^ be the set of all initial 

laws of the solutions belonging to lR[o,r] > which we regard as the set of initial distribution of 
mass that we are allowed to transport. 

Definition 2.9 (IR-MF). A Borel family {'n{^))xeR<i ^ {C{f),T];W^)) is said to be a 01- 
regular martingale flow (Ol-MF) (assoeiated to L) if the initial law of r]{x) is 5x, for every 
X G and, for every u G IRq; the probability measure f T/(x)dP(x) is a Oi-regular solution to 
the MP (assoeiated to L). 

We remark that we are not imposing that, for every x G Tf{x) is a IR-regular solution to 
the MP associated to T; the requirement is only in average, with respect to every admissible 
initial density h G IRq- Of course, from this condition and Proposition 2.8 we obtain that r]{x) 
is a solution of the MP, for u-a.e. x G for every u G IRq- For example, if we let lk[o,r] be 
the set of all narrowly continuous solutions of (3), then we operate no selection at all, and 
IR-MF’s are Borel selections {'n{x))x(zs,d of solutions of the MP, with r/(x) starting at 6x, for 
every x G The DiPerna-Lions theory is obtained if we let Ik be the set of all narrowly 
continuous solutions ut = of (3) with ||u||^^ < oo. 

We state (without implicitly assuming) some further properties of IR-regular solutions of 
MP’s and FPE’s that are useful in the next section. The first property is a stability property 
with respect to pointwise domination: for every F, i/, narrowly continuous solution of (3) such 
that, for some C > 0, 

h't < Cuf, for every t G [0,T] and iz G then D G (10) 

A useful property is stability with respect to convex combinations, i.e., for any i? G 

if, P-a.e. z £ Z, is IR-regular and (9) holds, then f ri,,dh{z) is Ik-regular. (11) 
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A reasonable converse should be stability with respect to disintegration, but there are several 
formulations: given any 3i-regular rj, writing for a regular conditional probability 

with respect to cq, we may require that 

for any i> G with z/ < Cijq for some C > 0, then J ri{x)iy{x) is 3i-regular, (12) 

or alternatively that 

for any i> £ JIq with i> <C r]Q, then J ri{x)iy{x) is 31-regular, (13) 

or even that 

for r/o-a.e. x G ri{x) is a fR-regular solution to the MP, (14) 

which is a rather strong condition: it formally implies the others whenever (11) holds true. 
Let us also notice that it does not hold when we deal with the DiPerna-Lions class introduced 
above, while (12) as well as (13) hold true. Moreover, an application of Theorem 2.5 shows 
that condition (12) is equivalent to (10). 

Due to technical reasons, we must introduce a slight extension of all the notions above, 
taking into account a family (3?[s^7’])5g[o,T]) where each consists of narrowly continuous 

solutions of the FPE associated to L, on [s,r]. Then, we let 3?^ be the set of all 1-marginals 
at time s for solutions belonging to 3l[^ and we refer to 3?-regular solutions of FPE’s and 
MP’s on [s,r], by natural extension of the definitions given on the interval [0,T]. We also 
assume that 


for any r, s G [0,r], with r < s, u £ then {i^t)t£[s,T] G (15) 

In particular, for any ly £ one has G 3?^. We also accordingly extend the notion of 

3?-MF by considering a family (^ 7 ( 5 , x))^g[Q where (^ 7 (s, is a fk-MF, for every 

s £ [0,T] (notice that we are not requiring joint measurability of (s,x) i-)- r]{s,x)). 

Remark 2.10 (Markov property). With the notation introduced above, we can state the 
Markov property via Chapman-Kolmogorov equations, for a tk-MF (??(s, x))^g[o,T],xeiR'^’ 

r]{s,x)t = J r]{s,y)t'r]{r,x)s, n-a.e. x £ , for every i? £ 3lr (16) 

for every r, s, t £ [0, T] with r < s < t. 

We obtain this property as a consequence of uniqueness, arguing e.g. as in [Fig08, Propo¬ 
sition 3.10]. However, let us remark that it could be be of independent interest to study regular 
flows that are also Markov, extending e.g. the approach in [SV06, Chapter 12]. Finally, much 
less is known about the strong Markov property for DiPerna-Lions flows, i.e., the validity of 
(16) with stopping times in place of deterministic times - perhaps one has to introduce some 
notion of “regular” stopping times. 

2.2 Equivalence betvizeen FPE’s, MP’s and flows 

The superposition principle provided by Theorem 2.5 allows for establishing a neat corre¬ 
spondence between “Eulerian” and “Lagrangian” descriptions, transferring well-posedness 
results both ways. Such a connection is firmly established in the deterministic case, see 
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e.g. [AC08, §4], and in the stochastic setting has been investigated e.g. in [Fig08, §2], in case 
of a DiPerna-Lions theory, and in [EK86, §4], for the classical theory (i.e., not in a selection 
sense). In this section, we provide a complete equivalence between well-posedness results for 
IR-regular solutions of FPE’s and MP’s. 

FPE’ s MP’s. Equivalence between existence result is straightforward, by lifting any 
solution V of the FPE, we obtain existence of solutions of the MP, so we focus on uniqueness. 
A simple result which transfers “uniqueness” is the following one: the non trivial implication 
ii) => i) follows from lifting two different solutions (see also [Fig08, Theorem 2.3]). 

Lemma 2.11. Let u G JIq. Then, the following conditions are equivalent: 

i) there exists at most one tH-regular solution v of (3) with uq = u. 

ii) ifr)^, rf are Li-regular solutions of the MP with ho = Vo = L, then Vt = Vt> t £ [0; T]. 

A stronger uniqueness result, for processes, can be obtained arguing as in [SV06, Theorem 
6.2.3] or [Pig08, Proposition 5.5]. Let us point out that here there appears a small gap with 
the deterministic literature, since a different argument [AC08, Theorem 9] shows uniqueness 
for MP’s assuming only (10), while we must consider also intermediate s G [0, T] (since the 
argument employed therein uses some conditioning which may not preserve the martingale 
property in general, but it does when the martingale is deterministic). 

Lemma 2.12 (transfer of uniqueness). Let IR = (iR[s,r])sG[o,T] satisfy (15) and (10), with s 
in place of 0, for s G [0,r]. Then, the following conditions are equivalent: 

i) for every s G [0,T] and u G there exists at most one u G iR[s,T] Vg = L. 

ii) for every s G [0,T], ifr]^, 77^ are Li-regular solutions of the MP on [s,T], with vl = Vs> 
then r}^ = rj^. 

Proof, ii) i). As in Lemma 2.11, n G with Vg = v we consider a (IR-regular) 

superposition solution ij- the uniqueness assumption entails that its 1-marginals are uniquely 
identified, i) ^ ii). The proof relies (implicitly) on the Markov property. Let s G [0,T] and 
T]^, rf be solutions of the MP on [s,T], with v) = vl- To deduce that = rf, we show that, 
for every n > 1 , the n-marginals of an rf coincide, i.e., for any s <ti < ... <tn <T and 
Ai, ..., An C Borel, it holds 

'rficti G Ai, ... ,et„ G An) = G Ai,..., G An). (17) 


We argue by induction on n > 1, the case n = 1 being a consequence of i) ^ ii) in 
Lemma 2.11 and property (15), i.e. we use the fact that {Vt)te[s,T] for i ^ {1)2} are IR-regular 
solutions, with vl = Vg- To perform the step from n to n -|- 1, we argue as follows. For fixed 
s <ti < ... < tn < tn+i < T and Ai,..., An, An+i C R*’* Borel sets, we let 


nr=iXT(eq) 

T 7 i(eii G Ai,... ,et„ G A„) 


:C'([s,r];R‘')^[0,oo), 


i.e., the density of conditioned with respect to fllLi ^ assume that the 

denominator above is not null: otherwise there is nothing to prove. Notice also that the 
inductive assumption gives (et„)j(yOT 7 ^) = {et„)fprf), since it amounts to (17) with A„ n B 
in place of A„, for every B C R'’* Borel. 
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For i G {1,2}, we let be the push-forward of the measure prj'^ with respect to the 
natural restriction from [s,T] to \tn,T], and notice that both are 3J-regular solutions of the 
MP on [tn-,T], with identical laws at 

by Lemma 2.6 and (15). By the implication i) ^ ii) in Lemma 2.11, we deduce in particular 
that thus 

±1± An+i) _ ^ +1 g ■^n+l) 

G An) £ Ai,... ,et„ £ An) 

hence we deduce the case n -|- 1 of (17). □ 

MP’s flows. In this case, both notions are “Lagrangian”, thus there is no need of the 
superposition principle here; most of the argument are just consequences of convexity and 
disintegration of measures. 

Although our actual well-posedness results are in the DiPerna-Lions case, where unique¬ 
ness is understood up to m-a.e. equivalence, where m is some “reference” cr-finite Borel measure 
on (i-e., m = Af^), for the sake of completeness, we provide a result assuming (14). 

Proposition 2.13. Consider the following conditions: 

i) for every n £ JIq, there exists a unique “Ji-regular solution rf to the MP, with initial law 
D, and the map i? ^ rj'^ is Borel; 

ii) for every i? £ IRq O'^d tR-MF’s {r]^{x))^^^d, {ri‘^{x))^^^d, one has r]^ = r/^, v-a.e. on 

Then, it always holds i) ^ ii), while ii) ^ i) holds true provided that some R-MF exists and 
both (11) and (14) hold. 

Proof, i) ^ ii) is straightforward, since regular conditional probabilities are essentially unique 
(a IR-MF is in particular a regular conditional probability of f rf(x)dP(x) with respect to cq). 

To show the implication ii) ^ i), let n £ Rq and ry^, ry^ be Ik-regular solutions of the 
MP, with initial law n. By disintegrating with respect to cq and (14) we may assume that 
n = 5x, for some x £ W^. Let {'n{x))x&R<^ be a Ik-MF (here we use the existence assumption) 
and define 

V'{x) ■■= X{xjtx}Vix) + X{x=x}V\ for i £ {1,2}, 
which are two different fk-MF’s since, for any p, £ fko, it holds, by (11), f T]^(x)dp(x) = 
f^x^x} V{x)dp{x) + p{x)r]^ £R. □ 

The result above is rather unsatisfactory in terms of existence of fk-MF’s, which seems a 
delicate problem, in general. For example, existence may follow if one assumes the validity 
of assumption i), (11), (14) and that Rq is a Borel of probability measures. Then, for every 
X £ such that G fko, there exists a unique Ik-regular solution of the MP, and by suitable 
definition for x not in such a set, we obtain a Ik-MF (which is then unique in the sense above). 
An easier existence result follows if we assume the domination condition 

for some cr-finite measure m, it holds P <C m, for every P G fko, (18) 

as in the DiPerna-Lions case. We also assume that m is minimal in the sense that, for every 
A C Borel with m(A) > 0, there exists some P G Iko with P(A) > 0. 
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Proposition 2.14. Let (18) hold, and consider the following conditions: 

i) for every i) G JIq, there exists a unique “Jl-regular solution rj'^ to the MP, with initial law 
V, and the map i? rj'^ is Borel; 

a) there exists a Ol-MF’s iT]{x))xGR‘^ IR-MF’s are m-a.e. unique, i.e. if {rj^ (x))and 
(j 7 ^(x))j,£Rd are Li-MF’s, then rj^ = rf , m-a.e. in M'^. 

If (11) and (12) holds, then i) ^ ii). If (11) and (13) are satisfied, then ii) ^ i). 

Proof, i) ^ ii). We have only to settle existence of some IR-MF, as uniqueness is trivial. For 
any probability F = nm G IRq we consider the unique IR-regular solution of the MP with 
initial law u and a regular conditional probability with respect to cq, {'n^{^))xGRd.- Then, for 
any um G it holds 

ri^{x) = r]^{x), m-a.e. x G X such that u{x) > 0 and v{x) > 0. (19) 

Indeed, it is sufficient to show that, for every e > 0 and every pm probability density, concen¬ 
trated on {rt > e, u > e}, with p uniformly bounded, it holds 


This, in turn, follows from uniqueness and (12); both members above are IR-regular solutions 
to the MP, with initial law pm < £~^Cum. 

Next, we notice that there must exists some wm G IRq equivalent to m, i.e., such that u> D 
m-a.e. in since m is equivalent to the supremum of all the measures in IRq (appropriately 
rescaled). Then, we define rj{x) := r]^{x), for x G M'^. To conclude that r]{x) is a 9J-MF, we 
use (19): given any probability um G h holds 

To prove ii) i), existence of 9J-regular solutions to the MP, given the existence of a 
9?-MF is trivial, so we focus on uniqueness. We let f], be a 9i-regular solution of the MP with 
some initial law and show that it must coincide with the one induced by the (unique) 9J-MF 
(j7(x))j,g]gd, i.e., f) = / rj{x)dfjo{x). To this aim, we let um G IRq be a probability measure 
equivalent to m, and consider the measure 


V 



r]{x)u{x)dm{x), 


which is a 9J-regular solution to the MP by (11), whose initial law is again equivalent to 
m. By disintegration with respect to cq, we obtain a Borel family of probability measures 
ivi^))xGR‘i^ which, by (13), provides a 31-MF and so by uniqueness it coincides with r]{x), for 
m-a.e. x G yielding 


r]{x) 


^dfjo{x) + ]^u{x)dm{x) 


= 7 ;V + 7 : ( r]{x)u{x)dm.{x) 


from which we conclude. 


□ 
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We end this section with some remarks on standard consequences of uniqueness: the 
Markov property and stability with respect to approximation. 

As in Remark 2.10, we consider IR-regular flows with respect to a family such 

that (15) holds. 

Proposition 2.15 (Markov property). Assume that uniqueness holds for “R-regular MP’s, in 
the sense that, for every s G [0,T], i? G fk*, there exists a unique “Jl-regular solution to the 
MP on [s,T], with initial law v. Then, for every tR-MF (^7(s, a:))sg[o,r],a;GR'*? (16) holds true, 
for every r, s, t G [0, T], with r < s <t. 

The proof is straightforward from the following identity between measures on ^([s, T]; M'^): 



ri{r, x)v{dx) 



j ri{r,x)sv{dx) 


{dy), 


which, in turn, holds true because both terms define fk-regular solutions of the MP on [s,T], 
with initial law f T](r, x)sh(dx}: this is obvious for the right hand side, while for the left hand 
side it is a consequence Proposition 2.6 and condition (15). 

Another well understood, but rather technical, property that sometimes follows from exis¬ 
tence and uniqueness is a non-quantitative version of stability with respect to approximations, 
which in this setting would read as follows. 


Proposition 2.16 (stability). For n > 1, let aP, fe” he Borel maps as in (2), let := 
T(a'^,6"') and let solve the MP associated to £"■. If 

i) there exists a unique R-regular solution rj of the MP associated to L = L{a,b) with 
Vo = h, 


a) it holds rjQ ^ n narrowly, a"' ^ a and ^ b pointwise as n ^ oo. 
Hi) for some convex, l.s.e functions ©i, 02 as in Theorem A.2 it holds 

Mr f-T r 


limsup / ei{\bf\)+ e 2 {\af\)dr]^dt < Qi {\bf\) + G 2 {\af\) dytdt, 


n—>-OD 


iv) every limit point in C'([0,T]; ^(M'^)) of {v^)n>i belongs to R, 
then r]^ ^ r] narrowly in 3^{C{[t),T];W^)). 

Notice that we do not require that r/"' are IR-regular: in general it does not even make 
sense, since IR is a class of solutions to the FPE associated to T, not to A proof of the 
result above would not be difficult, but it would require us to combine some technical results, 
such as those established in Section A.2 and [AC14, Lemma 23] to establish that (r/") is a 
tight sequence and any limit point provides a fk-regular solution to the MP associated to L; 
the conclusion is then straightforward from uniqueness. If (18) holds, then m-a.e. convergence 
in place of pointwise convergence of the coefficients is sufficient, if we also restrict to solutions 
whose marginals are absolutely continuous (as done, e.g. in [Fig08, Theorem 3.7]). 
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3 Well-posedness results 


In this section, we state and prove two results (Theorem 3.1 and Theorem 3.2) about exis¬ 
tence and uniqueness for solutions of the FPE (3), belonging to suitable classes of probability 
measures. In particular, as we are interested in the DiPerna-Lions theory, we deal with abso¬ 
lutely continuous with respect to the d-dimensional Lebesgue measure, satisfying 

some bounds on their density u : [0, T] x —)■ R. Besides such integrability conditions on 
the solution, we require (Sobolev) regularity assumptions on the coefficients a, b. 

In Section 3.1, we give some formal derivation of the energy estimates which eventually 
lead to well-posedness for FPE’s; in Section 3.2, we introduce the notation for Sobolev spaces 
and related basic facts; in Section 3.3 we state our results and compare their with some 
(related) existing literature; the technical heart of the matter is developed in Section 3.4, 
where crucial commutator inequalities are proved; in Section 3.5 we give a proof of main our 
results. 

3.1 Energy estimates and renormalized solutions 

As in the classical DiPerna-Lions theory (as well as in [EigOS]), we rely on energy inequalities 
satisfied by an absolutely continuous solution u = (ut)te[o,T] of (3) (i-e. Ht = Let us 

briefly sketch a formal derivation, where we assume all the quantities involved being smooth 
(solutions and coefficients). 

The main idea is to write the equation satisfied hy t i-> f /3(ut(x))dx, where /3 : R R 
is a smooth function (from a Lagrangian viewpoint, this amounts in choosing, as a test 
function, an expression involving the density of the solution u itself). The chain rule gives 
dtl3{u) = j3'{u)L*{u) and, by linearity, we consider separately the drift and diffusion terms. 
Straightforward calculus gives 

^'{u){b-V)*u = P'{u)V*{bu) = (3'iu)uV*b + P'{u)b-Vu = (6 • V)*,d(u) + [/3'(u)u -/3(u)] V*b. 

For the diffusion part, we first notice the identity (a : V^)*u = —V*(a • Vu) -|- V*((V*a)u), 
where V*a is the vector field (V*a)j = — Yl'j=i for i G {1,..., d}. Therefore, we obtain 

P'{u){a : = 

= - Y'iu)V*ia ■ Vu) + Y'{u)V*{{V*a)u) 

= - /3'(u)V*(a • Vu) + V*{{^*a)Y{u)) + [Y'{u)u - Y{u)] {V*Ya 
= - V*{a ■ VY{u)) - P”{u)a{Vu,Vu) + V*{{^*a)Y{u)) + [Y'{u)u - /3(u)] (V*)2a 
={a : VY*Y{u) - Y"{u)a{Vu,Vu) + [P'{u)u - Y{u)] {V*fa 

where (V*)^a = V*(V*a) = Y^ij dfjaij. Summing up, we have the identity 

dtYiu) = T*(/3(u)) - ^ Vu) - [/3'(u)u - /3(u)] (div T), (20) 

where divT := —V*6 — (V*)^a/2. By integrating over R'^, (formally T1 = 0), we deduce 

dt J I3{u)d^‘^ = - J ^-^aiVu, Vu)d^‘^ - j [/3'(u)u - /3(u)] (div T) 
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( 21 ) 


If 13 is convex with (3(0) = 0, so (3''{z) > 0 and I3'{z)z — I3(z) > 0, for z G M, we obtain 


dt 




(3'{u)u —/3(u) (div£) 


which is the key inequality we employ to show existence as well as uniqueness, for suitable 
choices of (3. For example, letting (3{z) = \z\~^, we deduce that if uq > 0, then ut > 0 for 
t G [0,T] (thus, for simplicity, we assume that Ut > 0 in what follows). In particular, to 
deduce uniqueness for solutions in L'(°(L((.) (for some r > 1), we show that the difference 
between any solutions u, v satisfies (21), with (3(z) = |z|^, and Gronwall lemma entails a 
uniform bound with respect to t G (0,T) for \\ut — Let us also notice that, with the 

choice (3(z) = |z|^, keeping track of the non-negative terms dropped above, we would obtain 
a bound for the “Sobolev energy” atiVut,'S/ut)d^^dt and, for /3(z) = \zY, with r > 2, of 
the energy 

r(r-l) r [ ul-^at(Vut,Vut)d^^dt = [ at(VuJ’^,VuT)d££^dt 

Jo Jwi {r/2 - 1)^ Jo J]Rd 

In the elliptic case, i.e. if there exists some constant A > 0 with a(v,v) > A |u|^ for every 
V G M'^, uniformly in (0, T) x we would deduce that any weak solution u actually belongs 
to the Sobolev space L^(IFi’^). Moreover, if we have no bounds on divL but only on (V*)^a, 
and b G L|(L“), we may still deduce some bound with respect to the energy z |zp, 

^ I \ut\^d^‘^, 

so that 

d^^ J \Vut\‘^d^^, 


dt I \ut\^d^^ < I \ut\^ ((V*)2at) 


+ 


4 life* 


A 


2 J utbt • Vutd^^ <^j 


and again Grownwall inequality leads to a bound for L^, uniform in t G [0,T]. Similarly, if 
r > 2, we use the inequality 2ab < + b'^ thus, for every e > 0, the term r f ul~^bt ■ Vutd^^ 

(assume for simplicity that u is non-negative) is estimated with 


r/2 - 1 


ud bt 




< 2 / iv«: 


r/2| 

t I 


'd^^ + 


\\bt 


2(r/2 - l)2e 


\ut\^d^^ 


and letting e = A, we may conclude again by a Grownall argument that 


sup \\ut\\Lr < 
iG[0,T] 


uoWlt exp 



LjL- 


+ 


2(r/2 - 1)2A 


\lIl- 


( 22 ) 


Let us finally remark that if we integrate (20) with respect to some function / G Al, with 
/ > 0 (instead of / = 1), we would deduce 


dt J f^(u)d^'^ < j(dt + Lf)(3{u)d^‘^ + J f [/3'{u)u - /3{u)] (divL)" d^'^. (23) 


The inequality above is so useful that weak solutions u of the FPE, which also satisfy 
(23) for every / G Al, / > 0, for (many) smooth convex functions (3, are called renormalized 
solutions [FigOS, Definition 4.9]. There are abstract results connecting well-posedness for 
FPE’s and the fact that every weak solution is renormalized, e.g. [FigOS, Lemma 4.10] (but 
see also [BC06] for a somewhat converse result, in the deterministic framework); here, for 
brevity, we limit ourselves to a direct proof of uniqueness of FPE’s from the validity of (20), 
e.g. with the special choice /3{z) = jz]”. 
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3.2 Sobolev spaces 

Before we state and prove our main results, we briefly introduce Sobolev spaces associated 
to the operators dt and £, together with some useful facts; we use throughout a compact 
notation extending that in Section 2.1. 

For Pi q ^ [l)Oo], the space is defined as the space of functions u G Lf(L^) such 

that the distributional derivative dtu is represented by a (unique) g G 

[ {dtf)utd^'^dt = - [ [ ftgtd^'^dti for every / G Ac- 

K'i Jo Jr'* 

We endow with the Banach norm + \\dtu\\]^P]^q^. A standard mollification 

argument, with respect to the variable t G (0,1), gives that A is dense in 1F/’^(L^), for p, 
q < oo (for a proof of this and the following results, we refer e.g. to [Sho97, §III.l]). In 
particular, the chain rule for dt extends to thus 

dtifPiu)) = {dtf)l3{u) + f/3'{u)dtu, for every f £ A, u £ 13 £ (M). 

Another straightforward consequence of the density of A is the fact that any u £ 
enjoys an absolutely continuous representative, i.e. there exists some u £ AC'^’([0, T]; 
such that ut = ut, for .if^-a.e. t £ (0,T). In particular the map: To(tt) := uq (trace at 0) 
is linear and continuous from wl’^{Lx) to L%. Moreover, t ut strongly differentiable at 
J^’^-a.e. t £ (0, T) and it holds -^u = dtu. 

We associate to the diffusion operator L some “Sobolev spaces”. An important role in 
our deductions is played by DP{L) (for p £ [l,oo)), defined as the abstract completion of A 
with respect to the norm ||/||£)p(£) := ll/UxiLg + which is well defined whenever 

a, b £ LjLx (actually, a more consistent notation for D'p{L) would be Ll{DP{Lt))). Let us 
remark however that, without further regularity assumptions, the extended operator D^{L) 3 
f I—7- Lf £ LI{Lx) may be multi-valued, but the assumptions on L that we impose in our 
results entail that the extension is single-valued. 

A useful fact is the following: if / G IF/’^(L§) n D^{L), then one can provide a sequence 
(/n)n>i ^ A converging towards / both in IF/’^(L§) and DP{L). Indeed, it is sufficient to 
consider first a sequence {gn)n>i ^ A converging towards / in DP{L), let p be a smooth 
probability density on M, and consider the approximation gn^m ■= dn * Pm (where we let 
Pm{t) = m~^p{t/m), f G M, and we carefully extend to a continuous function outside the 
set [0,T] X M'^). For (n,m) —)■ oo, the sequence gn,m converges towards / in DP{L)i because 
g >-3 g * Pm is a contraction in DP{L), as convolution with respect to t and the operator L 
commute; for fixed m > 1, the sequence gn,m converges towards / * pm, because g g* pm is 
continuous from into IF/’ (Lx), with norm smaller than ||pm|loo- Moreover, as m —)• oo, 

f * Pm converges towards / in wl’^{Lx), since / G IF/’^(L§) (this is exactly the standard 
mollihcation argument providing density of A in {Lx)). By a diagonal argument, we 
finally extract a sequence (/n)n>i as required. As a consequence, if n G (r > 1) is a 

narrowly continuous solution of (3), with a, 6 G LILx, then the weak formulation (5) extends 
to/GlF/’'(L”')nZI”'(£): 




[{dt + Lt)f] utd^<^dt = I fTUTd.^^^ - I huod.^\ 


(24) 
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where by fx G L^. and /o G L'^ we mean the continuous representative of / evaluated at T 
and 0 . 

Similarly, we introduce the space DP{L,aV 0 V) as the abstract completion of A with 
respect to the norm |||/| + \Lf\ + a(V/, V/)|| ri rP. Clearly, this is a space than DP{L), but 
is useful because the following chain rule holds, for D'P{L,aS/ (S> V), and 7 G C'^(M), with 7 ' 
and j" uniformly bounded: 

£(j(u)) = y(u)£(u) + {u)a{Vu, Vtt). 

As in the previous case, it might be that / £(/) and / 1 —)■ a(V/, V/) are not single¬ 

valued, but the identity above holds true with the natural interpretation (and in our results 
we introduce assumptions ensuring that these are well-defined functions). 

We always consider the divergence div£ be defined in the sense of distributions, i.e., as 
the linear operator AcB f I (provided that a, b are locally integrable). We 

say that divC G Ll{Lx) if there exists a (necessarily unique) g G Ll{Lx) such that 


'{o,r)xii 




'(0,T)xR'* 




for every f £ Ac- 


Similarly, divC“ G L\{L^) if, for some g G LI{Lx), the inequality < in place of equality above 
holds, for every f G A, with / > 0. If div£“ G L|(LS), then we can prove the following 
inequality, for u G D^'P{L,aV <81 V), and /3 G C'^(M) convex, with bounded derivatives as well 
as l3'{z)z — j3{z) bounded: 


[ Lil3'iu))ud^^+^ < [ [/3'iu)u-/3iu)]divL-d^^+^. (25) 

J{0,T)xM.<i J{0,T)xR<i 

Indeed, let p be a smooth convolution kernel on and consider the diffusion operator 
with smooth coefficients a * and b * p™ (where we let p™(x) = m~'^p{x / m)) . If we also 
assume u £ A, then the inequality above holds true by the derivation as in Section 3.1 above. 
The general case follows by approximation, letting hrst m —)■ 00 and then choosing vA £ A 
converging towards u m. u £ DP{L,a'V ( 8 * V). 

Besides these spaces associated with L, let us recall some features of standard Sobolev 
spaces and the smoothing properties of the standard heat semigroup (P“ )a>o on R'^. For 
p £ [ 1 , 00 ], we consider spaces 

IFi’^ := {f£Ll:Vf£ L^} , ^ ^ ^i,p . ^2 ^ ^ ^ 

endowed with the usual norms. 

A crucial fact for our deductions are quantitative inequalities for the smoothing effect of 
the heat semigroup (P“)a>o, which can be deduced by straightforward computations from 
the heat kernel in R*^. Of course, P“ is a contraction semigroup in Wx’^ as well as IFJ’^; 
moreover, integration by parts and Holder inequality give 

l|VP“/||iP < c||/||^P for every a G (0,-hoo), (26) 

with c depending on p G [l,oo] only (possibly also on the dimension d). Such inequalities, 
called — r in [AT14], play a fundamental role for our approach to continuity equations in 
metric measure spaces: their validity in abstract setups as well as in Riemannian manifolds 
follow e.g. from uniform lower bounds on the Ricci curvature. 
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Arguing similarly, it holds for p G [1, oo], i, j G {1,... d}, 

< c||/||£,p for every a G (0,oo). (27) 

Let us also notice that, as a 0, the left hand side in the two inequalities above are 
infinitesimal, for a standard density and uniform boundedness argument applies. 

Finally, another property that we occasionally use below is that, for p G (l,oo), one 
has Wx’^ := {f G Lx : A/ G Lx}, because of the Lx-boundedness of the second order Riesz 
transform / i—)■ V^A“^/, see e.g. [GTOl]. 

3.3 Well-posedness: statement of results 

We are in a position to state our main well-posedness results, which we split in two theorems: 
the first one deals with possibly degenerate diffusions, with Sobolev regular coefficients. 

Theorem 3.1 (degenerate case). Let p G (l,oo], r > 2p/{p— 1), and a, b be as in (2), with 

aGLliW^’P), bGLjiW^’P), and divL" G 

Then, for every probability density u G there exists a unique narrowly continuous solution 
u = (iit)tg[o,T] of the FPE (3) with uq = u and u G L^(L^). 

Actually, the technique employed provides (existence and) uniqueness even without the 
assumption that ii is a probability density. As a straightforward consequence of the result 
above and the equivalence established in the previous section, if we let IL be the class of 
narrowly continuous solutions u of the FPE (3), with u G Lf^(L^), we deduce existence and 
uniqueness for IL-regular martingale problems as well as for IL-regular martingale flows. The 
unique regular flow satisfies the Chapman-Kolmogorov equations (16). 

Our second statement deals with non-degenerate (elliptic) diffusions, i.e. if it holds, for 
some A > 0, a{v,v) > A|u|^, for every v G M^, a.e. in (0, T) x In such a case, we can 
remove one order of Sobolev regularity assumption from both coefficients, but we introduce 
Lipschitz regularity for 1 1 —)■ a*. 

Theorem 3.2 (bounded elliptic case). Let p G [2,oo], r > 2pf{p — 2) G [2,oo], and a, b be as 
in (2), with a G Lf'{L^) and elliptic, 

dtaGL^,,, aGL]{Wl^P), b G LlLf and {{V*fa)-GLl{L^). 

Then, for every probability density u G L^, there exists a unique narrowly continuous solution 
u = {ut)t£[o,T] of the FPE (3) with uq = u and u G Lf^(L^). 

Also in this case, as a straightforward consequence of the equivalence between Eulerian 
and Lagrangian descriptions, we obtain well-posedness for IR-regular martingale problems as 
well as IL-regular flows, with “Jl as in the previous case. 

Remark 3.3 (comparison with existing literature). The literature on the subject of Fokker- 
Planck equations and martingale problems is so vast and growing that we limit ourselves to 
a direct comparison only with very closely related and recent works. In particular, we stress 
some aspects which are different from the results appearing in [Fig08], [LBL08]. 
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In [LBL08], the approach is mostly Eulerian, dealing with FPE’s in divergence form 
dtut + V*{utb) = ^V*{aa*Vut), on (0, T) x 

with cj : ^ The main result in [LBL08] provides existence and uniqueness for the 

equation above, provided that 

beLliwlfl), V*6gL1(L-) 

To compare these assumptions, we must notice as in [LBL08, ^5.1] that with our notation 
a = aa* and the drift is actually b— ^V*a. In view of this correspondence, it might seem that 
Theorem 3.1 follows from their weaker assumptions: this follows in principle from a result 
of the type a := G if a ^ extending the well-known result [SV06, 

Lemma 3.2.3] that afl‘^ is Lipschitz whenever a G C'^. However, their conclusions are in 
fact weaker, and actually insufficient in order to obtain correspondent Lagrangian results: 
they prove existence and uniqueness in the class of narrowly continuous probability densities 
u G such that aVu G L^(L^); the latter (weak) regularity condition then prevents 

from a straightforward application of the results in Section 2.2. In conclusion, our result has 
(apparently) stronger regularity conditions on the coefficients, but draws stronger results and 
leads directly to well-posedness of regular martingale problems and flows. 

The problem arising from the condition aVu G L]{L‘],), which prevents a Lagrangian 
theory, is well understood in [Eig08], where much effort is put in showing, for the bounded 
elliptic case, uniqueness in the class of narrowly continuous probability densities u G 
[Eig08, Theorem 4-3]■ When compared with the assumptions of Theorem 3.2, an evident 
difference is that we require a first order condition a G Ll{Wx’^), while no such require¬ 
ment appear in [Fig08, Theorem 1.3], besides (with our notation) V*a, div£“ G Lff{L'ff). 
The technique we employ - approximation by the semigroup associated to the Dirichlet form 
f J a(V/, V/) - is the same as Figalli’s one, and in the elliptic case the novelty is more 
conceptual, providing a much cleaner derivation of commutator estimates, essentially by the 
same abstract arguments in the elliptic and the degenerate case. However, in the possibly 
degenerate case, our results are stronger, compare e.g. with ]Fig08, Theorem l.f], as we allow 
for much more general diffusion coefficients, and possibly unbounded terms - obtaining as well 
Lagrangian counterparts. 

In more recent years, further developments along these research lines appeared in the liter¬ 
ature, as well as different techniques (e.g. Crippa-DeLellis’ technique ]CDL08] was extended 
to SDE’s in fZhalO, RZIO]): of course, novelties and improvements appear in these devel¬ 
opments, but to the author’s knowledge that they address different aspects (such as strong 
solutions, equations with jumps, quantitative estimates, etc.) and there is no substantial 
overlap with our two results above. 

We also point out that the theory of measure-valued solutions (i.e., not necessarily abso¬ 
lutely continuous) Fokker-Planck equations, at least in the elliptic case, is well-developed and 
some results may be compared with ours. For example, (BDPRS07, Proposition 3.1] entail 
uniqueness if, for some p> d-\-2, a G LffHx^ is elliptic, b G L^Lx and 1at is Holder con¬ 
tinuous (locally uniformly in x). It is immediate to see that there is no inclusion between such 
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class of coefficients and that of Theorem 3.2, and in particular the hypothesis of our result are 
dimension-free (indeed, we are specializing a theory tailored for infinite dimensional spaces). 
However, the uniqueness class is smaller in our case, since we restrict from the very begin¬ 
ning to absolutely continuous solutions, which is nevertheless sufficient to entail a reasonable 
Lagrangian theory. Let us point out some recent developments [BDPR08, BRSll, BRS13] and 
in particular [BRS15] which also contains a survey of known results and methods for the de¬ 
generate case. Finally, we point out the monograph in preparation [BRKS15], which contains 
a detailed study and a vast bibliography on the subject. 

Let us briefly discuss some features of the two theorems above and their proof. First, exis¬ 
tence of weak solutions in the hypothesis stated above is a much easier task than uniqueness: 
for example, one can argue by approximation via convolution of the coefficients (and the ini¬ 
tial law) with a smooth kernel, so that the estimates on the coefficients are preserved, and one 
gains enough regularity (e.g. coefficients) so existence is available even at the Lagrangian 
level. Then, we have enough regularity so that the deductions which lead to inequality (20) 
apply and by a Gronwall argument we deduce a bound in in terms of divT“ only, 

and uniform in the approximation (in the elliptic case, we argue with (22) instead). By ex¬ 
tracting a weakly convergent sequence and by strong convergence of the approximations of 
the coefficients, we deduce that any weak limit point in L^(L^) is a weak solution to the FPE 
(3). In the elliptic case, we deduce as well existence for a solution u G L^(L^) n Lj (IFi’^). 
Let us also recall the approach [FigOS, Theorem 4.3], which is completely Eulerian (i.e., it 
relies on PDE’s techniques only), and has the advantage of yielding easily uniqueness, for 
solutions belonging to such a (smaller) space, which does not allow for applications of the 
theory developed in Section 2.2. 

In order to establish uniqueness of solutions, our aim is to rigorously establish (23) and 
(22), for (difference of) solutions u G As already remarked, the main problem is 

related to the regularity of u, in order to employ the standard calculus rules. Our strategy 
relies the well-known smoothing scheme, which dates back at least to [DL89]: for a G (0,1) 
we introduce some linear operator P“, acting on functions defined on (0,r) x such that, 
by defining := P“tt, we obtain an approximation of u sufficiently regular to rigorously 
obtain (23). Of course, the price that we pay is that u", in general, is not a solution of (3) 
and one has to carefully estimate the “error terms” thus appearing: our novel contribution 
indeed provides a systematic approach to such inequalities. 

To be more precise, in the cases that we consider, the operators (P“)q>o form a strongly- 
continuous Markov symmetric semigroup on L^((0,r) x so that, in particular, P“ 

preserves all L^Lx spaces, for p,q G [1, ooj. If we also prove that P" maps IF/’^(L”^) n D"^'(L) 
into itself, we may write, for / belonging to such space, 

r [ [{dt + Lt)f] ufd.^<^dt = [ ffiuTd^^- [ ffiuod.^^+ [ u [P“, (dt + Lt)] fd.^^+\ 
Jo J J J J{ 0 ,T)x'Rd 

(28) 

since the weak formulation (24) extends by density of A in W)’ (L’’ ) n (T). The commu¬ 
tator term appears as an algebraic way to highlight the identity as an equation for and 
all the issue is to show that it is infinitesimal as a 0. 

Next, we prove that P“ has a “smoothing effect”, in a sense that we can choose as 

a test function, and apply the chain rule with respect to dt and (25), so 

dt [ < [ [fi\uf)uf-ffiuf)]divLfid.^^+[ ut[P^,idt + Lt)]f3'K)d.^'^, 

J JiRd jRd 
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^^-a.e. t G (0,T) and in the sense of distributions on (0,T). Finally, we let a | 0, and by 
strong convergence of u" towards u in L|(L^), we are able to conclude, provided 


ut [P", {dt + Lt)] /3'(uf < e{a) ^ 0, 


in L^(0, T) as a I 0. 


3.4 Commutator inequalities 

In this section, we estimate the “error terms” involving the commutator between P“ and 
+ £. Our general strategy is a further development of that first introduced in [AT14], 
in the framework of continuity equations in metric measure spaces, and it is completely 
“coordinate free” and depends on an interpolation argument d la Bakry-Emery, namely 

J u [P“, dt + £] fd^'^ = j [P“u(at + L)f - u{dt + £)P“/] 

^ I lo dsd^^ 

= £ J [A, dt + £] P^-^fd^^ds, 

where we let A be the generator of (P“)o>o- It turns out that the commutator between A and 
+ £, reflecting the “relative regularity” between the chosen approximation and the target 
diffusion, depends upon natural quantities such as Sobolev regularity of the coefficients. 

In principle, this method provides very general results but, for the ease of exposition, we 
address separately only three cases, which are of particular interest: the case of a commutator 
between the Euclidean heat semigroup and a Sobolev derivation, which is a specialization 
of [AT14, Lemma 5.8] in the Euclidean case; that of a commutator between the Euclidean 
heat semigroup and a second-order Sobolev diffusion, which is apparently novel, that we settle 
by performing a “second order” interpolation argument; and finally that of the commutator 
between dt and a non-degenerate diffusion acting only the variable x G with t ^ at 
Lipschitz, which provides an alternative approach to Step 2 in [Eig08, Theorem 4.3]. 

We let throughout q G (1, oo], r, s £ (1, oo), with q~^ + r~^ -|- = 1 (one can deal with 

endpoint case at the price of more delicate approximations). 

Lemma 3.4. Let b G u G L^(L”) and f G It holds 



ut[P^M-^]ftd^^ 


dt < c||V 6 || 


L\L% 




LTLl 


LTL% ’ 


for a G (0,1), 


(29) 


where c G M is some constant (depending on the dimension d only). 


Actually, the proof below shows that V 6 can be replaced with the symmetric part of the 
derivative (also called deformation) := (V 6 -|- (V6)'”)/2, where r denotes the transpose 

operator. 

As a consequence of (29), the commutator operator L'^fWx’^) ^ f [P“, b ■ V]/ G L|L” 
extends to a linear continuous operator on Moreover, a standard density and uniform 

boundedness argument entails that, for / G 

[P “,6 • V]/ —> 0, strongly in Lj(L” ) as a | 0. 
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Proof. It is sufficient to argue assuming that b, u and / are sufficiently smooth, e.g., u G 
Ac., / G >1, as well as 6* G A, for i G {1,... ,d}, as the general inequality will follow by 
approximation (e.g. by convolution with a smooth kernel). Moreover, we argue at t G (0,T) 
hxed and then integrate over the interval (0,T): thus we omit to specify t G (0,T) in what 
follows. 

The curve s F{s) = f u^b ■ is then (7^(0, a), with 

^F{s) = j (Au")6 • V/“-" - u% ■ V(A/“-")dAf^ 

By straightforward integration by parts, we obtain the following alternative expression for 
the right hand side above: 


fL 

ds 


F{s) = / [((V6 + iyby)Vu% V/"-*) + (V*6) (Vn" • V/"-* - ii"A/“-")] 


If V*6 = 0, the conclusion is immediate, since we may estimate |T(q;) — ^(0)1 < Jg° 
and, by Holder inequality, 


/ 


u[P'^,b-V]fdFe^ 


< 2 


wy^b\y. Iivu^ 




Ils I|v/ 

r 2ds 


ds 
L% 


Jo \/s{a - s) 
<27r\\D^y-^b\yy\u\\^cJ\fhs. 

by (26) and using fj(s(l — s))~^^^ds = tt. 

The general case V*6 G is slightly more involved: let us hrst notice that the term 
(V*6)Vtt^ • V/"“^ can be estimated as above, adding a contribution vr||V*6||£^l to the in¬ 
equality. Finally, to estimate the contribution of (V*6)tt®A/““^ we do not put the absolute 
value inside integration with respect to s G (0, a), but exchange integration with respect to 
X and s, exploiting the identity 

r(v*6)n"A/“-"ds= -(v*b) ru^-^r-^ds. 

Jo Jo ds 

Next, to integrate by parts only “half of the derivative” with respect to s, we simply add 
{V*b)u°‘ times the quantity 


f-r- 


jL 

ds 


fas ^ 


thus 


ysdfc_sds 

ds 


< \u^{f-r)\+ r\{u^-uyAr-^\ds, 

Jo 


which, once integrated with respect to x G by Holder inequality and (27) is bounded from 
above by 

ds \ 


|V*6||i| \\uhr 


2 + c 


0 ys{s - a) 


This settles an analogue of (29) at fixed t G (0,T), and by integration with respect to 
t G (0,T), we obtain (29). □ 
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The constant c can be even independent of the dimension d of the underlying space, pro¬ 
vided that assume some bound directly on ||V*6||^I, and use a refined, dimension independent 
estimate for || : these are the key observation that lead to well-posedness on possibly 

infinite dimensional spaces, as developed in [AT 14], 

Lemma 3.5. Let a G Lj u G L“(L^) and f G L'^{Wx’^)- For a G (0,1), it holds 


ut[P^,at : V^]ftd^^ -a j ut[^,at : 


dt < c llV^al 


LfLl 11^ IILS°L 


(30) 


where c is some constant (depending on d only). Moreover, for u G L^(L^ n L|), it holds 


u[P“,a : V\P^u)d.^^ 


0, in L^{0,T), as a J, 0. 


(31) 


Proof. To establish (30), the underlying idea is to formally rewrite a : V^/ = a : (V^A“^)A/ 
and exploit the boundedness of the Riesz transform V^A“^ in L^, together with a second 
order interpolation along the heat semigroup. To make computations more transparent, we 
argue on coordinates, i.e., we fix i, j G {1,..., d} and consider the commutator 

[P“,a*’^a2^.]/ = - a*’^52/P“/). 

As in the proof of the previous lemma, we may also let u G Ac, f and a*’-^ be sufficiently 
regular, e.g. /, a*d G C^((0,T) x M*^), and argue at fixed t G (0,T). We consider the curve 


[0, a] 9 s !-)■ F(s) := 






which is C({0,a), with 


Ffs) 





u"[A,a*’^]52^/“-"dJ^^ 


since the Laplacian and partial derivatives commute. We write /i“ ^ := dFf^ ^ = {dFf)°‘ ^ 
(since derivatives and heat semigroup commute), let b := Va*’-^ and integrate by parts, ob¬ 
taining 







Differentiating once more, since F G C^{0,a), we obtain 


F"(s) = 2 


u*[A,6-V]/i“-"dA^‘^ + 


u"[A, (Aa*’^)]/r“-"dJ^'^. 


We introduce a second order interpolation based on the Taylor expansion 

rOL 

F(a) - F(0) - aF'(O) = / F"{a){a - a)da, 

Jo 

and we notice that the left hand side gives, up to integration on (0,r), the left hand side of 
(30). 
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Let us notice first how we would conclude in case V*6 = Aa*’-^ = 0, and then address the 
general case. As in the previous lemma, we obtain the identity 


and we estimate 


u"[A,6 • = -2 / ((v2a*’^')Vu", V/r“-") dAf" 


<^^=l|vV’^|L, 


V^s(a - s)3 


\LZ. 


Lt ■> 


where c is some constant. Integrating with respect to s G (0, a) and exploiting the factor 
(a — a) to compensate the bound the norm of we deduce (30). 

To address the general, we bound separately the terms 


j u"[A,6-V]/r° 


{a — s)ds and 


u"[A, (Aa*’^')]/i“-"dAf‘^ (a - s)ds. (32) 


To deal with former, we isolate a “leading term” which involves and we bound the 

remaining terms it by adding and subtracting suitable quantities, with the only difficulty that 
we must take into account the second order expansion. Precisely, after arguing as in the case 
Aa*’-^ = 0, we are left with estimating 


'0 




(33) 


and to this aim we add and subtract 


u“(Aa*’^)Ad“-"(a - s)dAf‘^ds = 


u“(Aa*’^)d?ijA^/“-"(a - s)d^^ds, (34) 


where we let Rijf '■= dfjA~^f be the second-order Riesz transform along the directions i, j. 
The difference between the (33) and (34) is easily bounded and to conclude, we exploit the 
identity 


/ A /“ ^(a - s)ds = / (a - s)a^/" ^ds =-/"aA/". 
do do 

and use the fact that the latter quantity is uniformly bounded (and that Rij is a bounded 
operator). 

To estimate the second expression in (32), we notice that 

thus we integrate by parts with respect to s G (0, a), 

£ £^{Aa^’^)dljr-^{a - s)ds = -au{Aa^’^)dl^r + £ u^{Aa^’^)dl^r-^. 

The first term in the right hand side above is bounded by c ||Aa*’-^ ||^5 II^IIl^ II/IIl®- write 

^ J u%Aa^'^)dljr-^d^‘^ds = J u^{Aa^’^)RijAr-^d^‘^ds. 
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To conclude, we argue once more by adding and subtracting 





u^{Aa^’^)Rijir - 


and estimating the differences involved. This settles the validity of (30), for smooth functions 
and at fixed t G (0,T). By integration and a density argument, the general case is deduced 
at once. 

Next, we prove (31), which follows from the fact that a f n[A, a : is infinites¬ 

imal, as a 4- 0: indeed, a standard uniform boundedness and density argument gives that the 
left hand side in (30) is infinitesimal as a 0. To show it, we initially argue in the case of 
smooth functions u, /, and for fixed i, j G {1,..., d}, we let b = and integrate by parts 




-2 j{b- Vu)dljrd^^ - j u{Aa^’^)dljrd^'^ 

-2 j {dljf)P‘^{b ■ Vu)d^‘^ - j u{Aa^’^)dljf^d^'^ 

-2 f (dljf) {[P“, 6 • V]u + (6 • Vu“)} d^^ - [ u(Aa*’^')A/“(iAf'^. 


Although the intermediate steps require some regularity for u, by the commutator estimate 
for Sobolev derivations established in the previous lemma, the resulting identity extends by 
continuity to u G L“(L^), / G Lf^{Wx’^)- Next, we specialize to the case / := tt“. By the 
strong convergence provided by Lemma 3.4 and uniform boundedness of adf jU°‘ in Lf°(L^), 
we have 


a 


j{dl^u^)[P^,bt^]utd^^ 


-G- 0, 


in L^(0, T) as a I 0. 


Similarly, it holds (recall that the left hand side in (27) is infinitesimal) 


a 


u(Aa*’^')A/“(iAf^ 


< ||Aa*’^| 


m 


LZ 


\\adl. 


.2a I 


U 


Finally, in order to handle the term a J(cl?jtt“)(6Vtt“)(iAf'^, the choice f = and the 
symmetry of a are crucial: we integrate by parts once, and since b = Va*’-^, we obtain 

f {dlju^)b ■ Vu^d^^ = - / [diu^{d],ka'’^)dku^ + (5iu“)(4a*’^')5|ju"] 

The first term, when multiplied by a, is clearly bounded and infinitesimal as a 0, so we 
focus on the last one. To show that it is bounded, we recall that a is symmetric and we are 
actually interested in bounds for the whole sum on i, j G {!,... d}; thus, by coupling the 
symmetric terms, it is sufficient to prove that 

is infinitesimal. This symmetric expression can be explicitly rewritten as 
1 1idka^’^)dk [(5,u“ + 9 ,u“)2 - (9,u“)2 - (5,u“)2] d^^, 


26 



and at this stage we integrate by parts once more, obtaining a bound in terms of 

which is sufficient to conclude (recall that the left hand side in (26) is infinitesimal). □ 


Finally, we deal with the bounded elliptic case; if a is bounded and elliptic, then the form 
L^{Wx’‘^) G / e-)■ Ja(V/,V/) is Dirichlet form, with associated Markov semigroup P" and 
(self-adjoint) generator A^/ = div(aV/), on its “abstract” domain D{Aa) (as given by the 
general theory of Dirichlet forms). When we choose P" as “smoothing operator”, the main 
difficulty is to prove that it preserves regularity with respect to t G (0, T), thus we need some 
estimate for the commutator [P“,i9t], which we initially define in following the weak sense, 
for u G ylc, / G yi: 


'(0,T)xI 




'{ 0 ,T)xR‘^ 


[{py)dtf + {dtu)P2f]d^^+<^. 


Lemma 3.6. Let a be bounded and elliptic, with dta G L“(L“). Then, for every a G (0,1), 
u G Ac, f G A, it holds 


J uK,dt]fd^^+'^ 


U 


LlLl 


LlLl > 


(35) 


where c is a constant (depending only on the ellipticity constant X). 

Thanks to this lemma and a density argument, for / G fF/’^(L^), we deduce that P"/ G 
VF/’^(L^), and the “strong” commutator [P",9t]/ := Pfdtf — dtPaf is well defined and it 
belongs to L^(L^). Moreover, the usual uniform boundedness arguments shows that, for 
u G Lj(L^) and any family {fa)a>o T L‘({Wx’‘^) converging in it holds 


—y 0, SiS CK 0. 

7(o,r)xR^ 

Proof. We provide the following analogue of (35), where dt is replaced by — I), where 

T'^/(t, x) = f{t + cr, x), and I is the identity operator (we also choose ct / 0 small enough, to 
avoid boundary terms, thanks to the assumption u G Ac)' 


'(0,T)xI 




< c 


L?°L° 


U 


L^tLl 




(which c depending on A only). Once this is is settled, we may let ct —)■ 0 and pass to the limit 
in the weak formulation. Let us notice that the identity operator plays no role above, and 
everything reduces to estimate J u[P“, By first-order interpolation along the 

semigroup P*, for s G (0,q;), it is sufficient to bound the infinitesimal commutator 

J u^[Aa,T^]r-"d.^^+^ = J ((T'"a)Vu",VT'"/“-")) - (aVu^VT'"/“-") 

where we performed integration by parts with respect to the variable x G and the change 
of variables t t + a in the first integral. We have therefore the bound (we are actually 
interpolating also along the semigroup u —)■ T'^) 





J (ra)(Vu^VT'^/“-")(i^^+'=' 


dr 
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which by drT^a = T^dta gives the thesis, after an application of Holder inequality and using 
the smoothing effect in Lf(L^) of Pa, i.e. ||Vu^||^ 22^2 < II^IIl 2 j;^ 2 . □ 

It would be natural to extend the argument above for more general exponents beyond the 
case above; the main issue being that a smoothing effect for Pa akin to (26) is not ensured 
by Dirichlet form theory, when the exponent involved is different from 2. It seems plausible 
however to replace with and require only dta G L\{L'^) (as the semigroup 

acts only fiberwise). 

Remark 3.7 (trace semigroup at f = 0). Another consequence of Sobolev regularity of the 
lemma above is existence of a “trace” semigroup, e.g. att = 0, defined as follows: for f G L^, 
consider a constant extension f{t,x) = f{x) for {t,x) G (0,T) x and let Pg/ be the trace 
of the Sobolev function P"/ at t = 0. Alternatively, this can be obtained as the semigroup 
generated by the bilinear form given by the trace at 0 of a. 

3.5 Proof of well-posedness results 

In this section, we address the proof of Theorem 3.1 and Theorem 3.2. As already remarked, 
existence is easily settled by approximations, so we focus on uniqueness. 

Proof of Theorem 3.1. Let u be the difference between any two narrowly continuous solutions 
in L“(L^) and let P" be the heat semigroup on R'^, extended on (0, T) x R^ by acting on each 
fiber {t} x R*^, t G [0,T]. For a > 0, P" maps IT/’^(L^') n D'^'{L) into itself, as /" := P"/ 
is C| with respect to the variable x G R'^, for .if^-a.e. t G (0, T) (to approximate P"/ with 
functions in A, we argue by convolution with a smooth kernel with respect to f G (0, T)), 
thus (28) holds true for / in such a space: 

r [ [{dt + Lt)f] ufd^<^dt = [ f^UTd.^^- [ f^uod^'^A [ u [P“, {dt + Ti)] fd.^^+'^. 
Jo J J J J(0,T)xR‘i 

(36) 

For a > 0, we also have u°‘ G D” (T): to use u" as a test function, we deduce that 
u°‘ G which follows directly from the equation satisfied by Indeed, (36) for 

f G Ac entails that the distributional derivative dtu^ coincides with the distribution L*u°‘, 
which is represented by a function, namely 

= V*(6u") + i(V*)2(au“) = -(divT)u“ + Lu^ + ^(V*a) • (Vu“) G 

Therefore, cltu“ G Ll{L^), and admits an absolutely continuous continuous representative, 
which must coincide with the one that we would obtain by acting directly to the narrowly 
continuous representative ut with the heat semigroup P“, at every t G [0,T]: it holds in par¬ 
ticular Uq = 0, since ug = 0. Moreover, the curve t e-)- J^d{uf )‘^dA^’'^ is absolutely continuous, 
with distributional and .if^-a.e. derivative ^ f = 2 J {dtu°')u°‘d.^*^. 

We are in a position to let u°‘ in the weak formulation (36), to obtain 

r [ [(5i + Lt)u^] ufdAT'^dt = [ + [ u [P“, Lt] 

Jo J J J{o,T)xmd 


28 



If we choose instead a test function t i—)■ f{t)uf, with / G C^[0,T) and we apply (25), we 
eventually deduce the inequality 


4: [ < ||div£,-|| _ [ + [ ut [P“, (dt + £i)] 

at J J^d J^d 

^^-a.e. t G (0,T) and in the sense of distributions on (0,T). Gronwall lemma gives 
||it"|lioci2 < exp|||div£"||^i^,„| / [ utiP"^, {dt + Lt)]ufd^‘^ dt. 

As a consequence of Lemma 3.5, we deduce ||rt||i~ 2.2 ^ Uminfcj^o =0. □ 

Proof of Theorem 3.2. In our smoothing scheme, we choose P" = P“ be the semigroup asso¬ 
ciated to the Dirichlet form f J a(V/, as introduced in the previous section. 

A first step consists in showing that (36) holds true, and we see it as a consequence of the 
fact that P" maps IL/’^(L^) r\D'^{L) into itself; if / G VL/’^(L^), then Lemma 3.6 shows that 
/“ G iy/’^(L^) as well; to show /" G D'^{L), we rely on the assumption on a G Lj{Wx’^), 
and show that the smooth approximations obtained by means of the standard heat semigroup 
P^(/“) converge towards /" in D‘^{L), i.e. £P^(/“) —)■ £(/") in Lj{L‘^) (this is the only point 
where we use the first order regularity assumption on a). Such convergence can be seen by the 
commutator lemma for Sobolev vector fields. Lemma 3.4, noticing that the claim convergence 
amounts to show 

[L,Pir ^Oin LliLl), 

but since derivatives and the standard heat semigroup commute, it holds 

d d 

[£, pir = E ^ 0 

i,j=l i=l 

since djf°‘ G and Lemma 3.4 shows convergence towards 0 in LIL^, as s 0. 

As a second step, we notice that we may let rt" be a test function in (36) indeed, it holds 
ii" G H^’‘^{L, a{V (8> V)) by what we just proved, while the fact that dfU^ is represented by 
some function in follows from a duality argument: for a.e. t G (0, T) the linear functional 
/ J^d Utt^tf^ is bounded in L^. From (36), we have 

j [iuf)\V*)^at + utbtVuf + ut K,dt] uf]d.^‘^, 

where we applied (25) only for the diffusion part a : V^, as we deal with the drift term 
separately, using the inequality 

\utbt ■ Vu“| < A |Vn“p + 4A-1 \ut\^\bt\^ , 

to bound the contribution of the drift part. To conclude, we apply Gronwall inequality and 
finally let a j, 0, using (31) to deduce that the commutator term gives no contribution in the 
limit and uniqueness holds. □ 
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A The superposition principle for multidimensional diffusions 


To prove Theorem 2.5, we follow a general scheme, whose structure is shared by many proofs of 
superposition principles appearing in the literature, see e.g. [AGS08, Theorem 8.2.1], [AC08, 
Theorem 12], [AF09, Theorem 4.5], [Fig08, Theorem 2.6], [AT14, Theorem 7.1], that we sum¬ 
marize below. The derivation is rather elementary, although the “right” underlying framework 
would that of Young (or random) measures. For simplicity, we let T = 1 in this section. Let 
^ ^ be a narrowly continuous weak solution of the FPE (3). To deduce 

existence of a superposition solution for z/, we perform the following steps. 

Step 1 (approximation). We build from v a sequence of solutions (z^"')n of FPE’s associated 
to diffusion operators (£"■)„, for which the superposition principle is already known to hold, 
thus obtaining a sequence of superposition solutions (r/"')n of MP’s. Here, the difficulty is to 
exhibit a sufficiently good approximation, so that v'^ converge towards zz, e.g., narrowly, and 
T” towards T, in a sense to be made precise, as n —)• oo. 

Step 2 (tightness). We prove that (r7")n C .^(C'([0,1];M'^)) is tight, yielding a narrow limit 
point rf. By Ascoli-Arzela criterion, this step reduces to show uniform bounds on the modulus 
of continuity of the canonical process (et)ie[o,i] with respect to r/”. 

Step 3 (limit). From convergence zv” — )■ u, ^ T, as n — )• oo, we conclude that rj is 
a superposition solution for u. Here, the problem is to deal with convergence for possibly 
non-continuous functions, as they involve the coefficients a, b. 


A.l Approximation 

We approximate the limit solution by means of mollification by convolutions or push-forwards 
via smooth maps (in probabilistic jargon, by conditioning with respect to some observables). 

Push forward via smooth maps. This technique is inspired by the approach in [AT14, 
Theorem 7.1]. Let tt G C'^(M‘^;M'^) TT = (vr^,..., tt'^), with uniformly bounded first and 
second derivatives. Then, it is possible to define a diffusion operator 7r(T) on such that 
TTjjZ/ := (vrjjZ/t)j£[o_i] is a solution to the associated FPE (in duality with A = C'j|’^((0,T) x R'^)). 
Indeed, the composition / o 7r(t, x) := f{t, 7r{x)) belongs to A, and if we let / o vr in the weak 
formulation (5), the chain rule gives 


L{f OTT 


d 1 ^ 

) = ^T(7r*)[(5i/)o7r] + - ^ a(V7r*, Vtt-^) [(9ij/) o tt] . 
2=1 ^J = l 


We define, for (t, x) G (0, T) x R'^, 


7r{a)l’^(x) := a(V7r*, Vvr-^) | tt = x 


dvTfj a(V7r*, V7r'^)z/t] 
dzrpt 




for i,j G {1,.. .,d}, 


TT(b)l{x) := Et-t [^(vr*) | vr = x] = ——^-^(x), for z G {1,..., d}. 

dn^ut 

Then, 7r(T) := £(7r(a),7r(b)) is a diffusion operator on R*^ and the (narrowly continuous) curve 
of measures TVjjiy is a weak solution of the FPE dt7r(iy) = tt(L)*tt(v), in (0,r) x R'^. Let us 
remark that 7r(T) depends upon z/, although it is not evident in the notation. 
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Since conditional expectations reduce norms and the derivatives of tt* are uniformly 
bounded, integral bounds on a, b are naturally transferred on 7r(a), 7r(6) (in particular, uni¬ 
form bounds). However, local integrability conditions could be not preserved. 

Mollification by convolutions. This is a more standard technique, already employed 
e.g. in [AGS08, Theorem 8.2.1] and [Fig08, Theorem 2.6]. Let p > 0 be a smooth probability 
density (with respect to.if'^), with full support. Then, the family of measures u * p := 
* /o)te[o,i]) solve a FPE associated to a suitably defined diffusion operator. Indeed, for 
f G A, it holds f * p G A with 


d 1 ^ ^ 1 ^ 

^if*P) = Yl b"di{f * P) + X X] * P) = X] b'idif) *P + i^Y^ * P 


i=l ^=1 

since derivatives and convolution commute. We define 

d{A'^vt) * p _ d{Uvt) * P 




: = 




, for i,j G { 1 ,... , 4 . 


d{vt * p) ’ ' d{vt * p) 

so {ut * /o)tG[o,i] is a weak solution of the FPE associated to := T(a^, 6 ^), as 


dt 


fd{iy * p) 


j {dtf)*pdv = 


dt{f * p) du 


T(/*p) du 


J LPfd{u*p). 


Integrability and regularity properties of and are collected by the following lemma, 
see [AGS08, Lemma 8.1.10] for a detailed proof. 

Lemma A.l. Let p be a smooth probability kernel on with p > 0 everywhere and |VV| < 
Cp, for i G {1,... k}, for some eonstant C > Q. Let p, u G with p u. Then, it 

holds p * p u * p, and the funetion 

d(u *p) J p(x - y) duly) 


provides a C^{W^) version of the density d{p* p)/d{u * p). Moreover, for every eonvex, lower 
semieontinuous function 0 : M ^ [ 0 ,00], it holds 


/ 


0 


f d{p*p) \ 
\d{u*p)) 


d{u * p) < 



(37) 


Similar conclusions hold when p = {pt)tG[o,i] ^ is a Borel curve and u = 

(^i)iG[o,i] ^ is narrowly continuous, with pt <C ut for every t G [0,1]. In addition, it 

holds 


sup 

te[o,i] 


djpt * p) 
d{ut * p) 


Ct(B) 


< 00, 


for every open bounded set H C R'^. 


When applied to a solution u = (t't)tg[o,i] of the FPE (3) we deduce that, if a, 6 G L'p{u), 
then aP, bP G LP{u * p) (for p G [l,oo]) and a^, b^ are C'^(R‘^), uniformly in t G [0,1], with 
uniformly bounded first and second (spatial) derivatives on compact sets of [0,T] x R*^. 
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A.2 Tightness 

We prove a compactness criterion for solutions of martingale problems, under minimal in- 
tegrability conditions on the coefficients. In the deterministic case, tightness is achieved by 
estimating the metric velocity of absolutely continuous curves which solve the ODE; in the 
stochastic case, we rely on analogous results for martingales, using Burkholder-Davis-Gundy 
inequalities and an argument reminiscent of Levy’s modulus of continuity for the Brownian 
motion, to estimate the modulus of continuity of the canonical process (yielding in some cases 
Holder regularity). 

Theorem A.2. Let 6, ©i, ©2 : [0,+oo) —)■ [0,+oo) he functions with ©i, ©2 convex, l.s.c., 


lim 9{x) = lim 

X—¥00 X^ + OO X 


lim 




©2(3^) 

X 


= 00 


and, for some constant C > 0, © 2 ( 20 :) < C'© 2 (x), for x > 0. Then, there exists some 
coercive function 'L : (^([ 0 , 1 ];M) —)■ [ 0 ,+oo] such that, for every filtered probability space 
(D, (3't)^g[oq], P) and progressively measurable processes ip = {ift), fi = « = {at)t with 


[0,1] 9 t !-)■ Mt := pt — [ fis ds, and [0,1] 9 t e-)- — 

Jo 

P-a.s. continuous local martingales, and a > 0 P-a.s. it holds 


a, ds 


< E 


9{po) + [©i(|/3i|) + ©2(at)]dt 


(38) 


Proof. We prove separately the existence of functions 'I'l and ^ 2 ; depending respectively on 
© 1 , ©2 only, taking integer and non-negative values, such that, if we let 


^(7) ;= 6I(|7o|) inf {^1(7^)^'2(7^)} , ( 39 ) 

7^-1-7^=7 


then 'k is coercive and (38) holds. For every e > Q, i £ {1,2} we let 5i = 5i^e be the largest 
number in the form 1 /n, (with n > 1 natural) such that Qfie''/5)5 > such a choice is 
possible because lima,^+oo &i{x)/x = 00 . Then, we introduce the closed sets 

A(e) := (7 e <^([0,1],I^) : sup sup |7s - 7(fc-i)5j < e|, (40) 

^ k=l,...,S~^ sG[(fc-l)(5i,fc(5i] ^ 

and we let ^*( 7 ) := Em>o("^ + ^)XAi (2 — m'^c . By construction, ^( 7 ) < m entails 7 G A(2 ^), 
for every k > m. 

To show that T defined by (39) is coercive, it is sufficient to apply Ascoli-Arzela criterion, 
noticing that 7 G {T < mj can be decomposed as the sum of two curves 71 + 72 , and 7 * 
(i G {1,2}) admits the following modulus of continuity 


J 2^ ^ if X G A: > m, 

{ if X G +00). 


(41) 


To show that (38) holds, we assume that the right hand side therein is finite. The as¬ 
sumptions entail therefore that {Mt)t is a P-a.s. continuous local martingale, whose quadratic 


32 



variation process is t J^agcls. If we let 7 ^ := Jq Psds and 7 I = Mt, for t G [0,1], then the 
left hand side in (38) is smaller than 

® + ^ 1 ( 7 ^) + '^' 2 ( 7 ^)] ^ IE [^(y’o)] + {fn + 1)E \_Xai(2-'^) 07 ^+ Xyi2(2-’") o 7^] • 

m>0 


Next, we focus on the addends in the series above, writing for brevity e in place of 2 ™'. For 
i G {1,2}, using (40), we have 


^[XAi(eY°l'']= p( sup 


(7')E > e) 


o. 

<E^((7*)fc>e), 


where we write, {y% ■= snp,^i(^k-i)Si,kSi] H “ 7(fc_i)<5j- 

Let us focus on the case i = 1 (thus we write <5 = di, 0 = 0i). Since | 7 s — 7 t| < f* |/?r| dr, 
we estimate 


^(( 7 ') 


E 

® (7 I(k-i)s 1/5^1 ds) 


pkS 

>e)<- 

e{s/5) 

- < eE 

/ 0 (1/3.1) ds 

J(fc-1)5 


where the last inequality is a consequence of Jensen’s inequality and our preliminary choice 
for J. Summing upon k G {l,..., }, we conclude that E ['L o 7] < c E [0 (|/ 3 <i|)] ds, for 

some constant c > 0 (in this case, the constant does not even depend upon 0). 

To deal with the case i = 2 (again, we omit to specify i in what follows), i.e., the martingale 
part, for each /c G {l,..., we estimate from above, 


P{Mt>e)< 


E[0((M*)VJ)] 

0(e2/(5) 


E 


< ce- 


0 


1 rkS 

5 J(k-l)S 


a.ds 


0 (e 2 /J) 


where c© is some constant depending on 0 only: indeed, it is sufficient to apply Burkholder- 
Davis-Gundy inequalities, e.g. in the form [LLP80, Theorem 2.1], to the martingale Mg := 
s G [0,(5] and the convex function with “moderate growth” x i-G 0(x2). By 
Jensen’s inequality and our definition of Sg we conclude that 


E 


0 


/ 1 rkS 

1^5 J(k- 


1)5 


a.ds 


0 (e 2 /J) 


< eE 


rkS 


'(k-l)S 


0 (q;^) ds 


As in the previous case, by summing upon fe G {O,..., (5 ^}, we deduce that 


IE [xA(eY o M] < ec©E 


Uo 


0 {as) ds 


and so we deduce the desired bound for E [4'(M)]. 


□ 


Corollary A.3. In the situation of the theorem above, let 0i(x) = and 02(a:) = 
for Pi, p2 G ( 1 , 00 ). Then, for every r > 0 with r < r{pi,p2) '■= min|l — \ 

holds 

P I limsup sup E—E = 0 ) =1. 

y HO |i-s|<h 7-^1 J 
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Proof. It is sufficient to 5i := for i G {1,2}. Thanks to this choice, the probabilities 
of Ai{2~^Y decay sufficiently fast as m —)• oo so that, by Borel-Cantelli lemma, there exists 
P-a.s. some m > 1 such that the curve (v?t)tG[o,i] can be written as a sum of two curves 
having dehned in (41), as a modulus of continuity. This entails r-Holder estimates for 
(p: since the condition on r is open-ended, thus arguing with a f slightly larger than r, the 
thesis follows. □ 


It is not clear if p in the previous lemma is actually P-a.s. Holder continuous with exponent 
r{pi,p 2 ). one might exploit the existence of functions 0* (i G {1,2}) with 0j(x)/ —)• oo 

as X —7- oo, and the right hand side in (38) still finite, but it seems not sufficient. 

Corollary A.4. Let a, b be Borel maps as in (2), let rj G .^((^([0, 1]; M'^)) be a solution of 
the martingale problem associated to L{a,b). For any 6, 0i and 02, as in the theorem above, 
let T be the associated coercive functional. Then, for every f G C'j|’^((0,T) x R'^), it holds 


'^{ftoet)dr]< / e{\fo\)dr]o+ / [0i(|Tt/|) + 02(at(V/t, V/t))] dr/tdf. 


Proof. We prove that t pt ■= ft ° satishes the assumptions of Theorem A.2, with 
(3t := {dt + Lt)ft o et and at := at(Vft, Vft) o et. For simplicity of notation, we omit to write 
et below (since its appearance is quite natural). Since both / and G C'^’^((0,T) x R'^), 


both 


t mI := ft- fo- I {ds Ls)f ds, and t M/' ■.= f'f - f'f - I {ds + Ls)f ds 


r ._ f2 


are martingales. By developing we see that 

t e-)• (m/~ y 4" ds — 2 J{9s + '^s)f 


{dr Lr)fdr - ft 


ds 


is also a martingale. We add and subtract the process t e-)- 2 fs{dt + L)fsds, thus 


ft ft 

tr^{M/)‘^— / asds + 2 {ds + Ls)f 

Jo Jo 


ft- fs- / {dr + J^r)fdr 


ds 


is a martingale. To conclude, we notice that 


11 —7- / {ds + Ls)f 


ft - fs - / {dr + Lr)f dr 


ds 




ds 


is a local martingale (see also [SV06, Theorem 1.2.8]). Indeed, by a stopping time argument, 
we are easily reduced to the case where is replaced by a martingale M with Mi G L°°(P), 
thus /q Ps {Mt — Mg) ds G T^(P), for t G [0,1]. To prove that increments are orthogonal, we 
let t G [0,1] and show that 


E 



/3s {Ml - Ms) ds I Tt 



Mg) ds. 


By the integrability assumptions, we exchange between conditional expectation and integra¬ 
tion. The thesis follows by direct consideration of the cases, s G [0,t] and s G (t, Ij. □ 


34 



A.3 Limit 


In the third step, we assume that the probability measures obtained as superposition 

solutions for a suitable approximating sequence narrowly converge in 

towards some limit rj. The fact that ij provides a superposition solution for i/ is not straight¬ 
forward, since we must deal with a limit in the weak formulation, where terms involving the 
coefficients a, b appear (in general, they are not continuous). 

Indeed, r] G ^((^([0,1];M'^)) is a solution of the martingale problem associated to L{a,b) 
if and only if the following property holds: for every s, t G [ 0 , 1 ] with s < t, for every 
/ G 1] X (with 11 / 11 (^ 1,2 < 1) and for every bounded continuous and T^-measurable 

function g on (^([O,T];(with Us'Hoo < 1) h holds 


ftoet- fsOCs- / [{dt Lr)f] o Cr dr 

J S 


As the correspondent identity holds for rj^ and T”, i.e. 


J 9 ft o et - fs o es - [{dt + T”)/] o dr 


dr] = 0 . 


drj^ = 0 , 


to deduce that ij is a solution to the martingale problem associated to L, since / and dtf are 
bounded and continuous, the crucial limit is 




(£”/) o Sr dr 


dr- 9 


{Lrf) o Br dr 


dr] —)■ 0 , 


(42) 


whose validity we now investigate, according to the approximations from Section A.l. 

Push forward via smooth maps. For n > 1, let tt” G C^{W^,W^) with {TTn)n converg¬ 
ing to the identity map locally uniformly and assume that the sequence of first and second 
derivatives converge (towards the respective limits), pointwise and uniformly bounded, i.e., 
V 7 r"'(x) — 7 - Id for x G V^ 7 r"'(x) —)■ 0, for every x G and there exists some constant 
C > 0 such that |V* 7 r"'(x)| < C, for x G and i G {1,2}. 

Let z/” = 7r"'(£), and let r]^ be corresponding superposition solution. To prove that 

any narrow limit point r] is indeed a superposition solution for z/, with respect to the diffusion 
operator L, we add and subtract the term 


(Lrf) o Br dr 


dr- 9 


irf) 


o e, 


dr 


dr] 


in (42), where L = L {a, b) is any diffusion operator on whose coefficients a, b are contin¬ 
uous and compactly supported. The difference terms above are infinitesimal as n —)• oo, by 
narrow convergence of 77 "', thus we estimate (42), as n —)• 00 , in terms of 

limsup f \L"'f — Lf\di^^ + ( \Lf — Lf\ d7T{iy). (43) 

n^oo J J 

Let us focus on first term above, at fixed n > 1 (for simplicity of notation, we drop the 
dependence upon n). Recalling the definition of 7 r(T),integration with respect to the push- 
forward measure gives 

J K(^)/ - ^/| = J [^(/ o vr) I tt] - (£/) o 7 r| du. 
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Being (Lf) o n a function of tt, up to i/-negligible sets, we have 

j [£(/ o vr) I tt] - (Lf) oTr\du = j \Ey [L{f o tt) - {Lf) o tt | tt] | di/ 

< j |^(/o7r) - {Lf)oTr\du, 

since conditional expectation reduces L^(i/)-norms. Writing explicitly the difference 
-• k k 

^(/o7r)-(£/)o7r = 2 X] [«(Vvr*, Vtt-^) - o tt] {di,jf)oTT+ -V ott {dif)oTr, 


i,j=i 


i=l 


and recalling that ||/||(^i,2 < 1, we conclude that 


TT{L)f - Lf \ d-K^v < / X ^ |«(V7r*, Vvr-^) — a*’-^ oTT\dv + E \L{^^)-h 


— 6* O TT 


*J = 1 


i=l 


dv. 


Letting n —)• oo (recall that tt = tt” above), using the assumption on the convergence of 
towards the identity map (in particular, we use Lebesgue dominated convergence w.r.t. the 
measure v), we deduce that (43) is bounded from above by twice 


k k 

/ ^ ^ 1“*’'^ “ a*’^| di/ + J '^\b" - V\ du. 

i,j=l 2=1 


To conclude, we choose d, b that minimize the right hand side above: this can be made 
arbitrary small, by the density of continuous and compactly supported functions in L^(z/)). 

Mollification by convolution. In this case, the argument is similar, and it is more 
standard, see e.g. [AGS08, Theorem 8.2.1], thus we only sketch it. Given a sequence p"' of 
probability densities on such that —)■ do narrowly as n —)• oo, let zz” = zz * p"' and T” 

be the diffusion operator introduced in Section A.l. We add and subtract, in (42), 


Lrf O Br dr 


r ft 


dv^- 9 


Lrf o Br dr 


dr], 


where L = L{a, b) has continuous and compactly supported coefficients. Let tJ be a common 
(bounded and continuous) modulus of continuity for d, h. 

As in the previous case, narrow convergence implies that the absolute value of (42) is 
bounded from above, as n —)• oo, by 


limsup f \L^f — Lf\di'^ + 
n^oo J 


\Lf-Lf\du. 


First, we prove that lim„_j.oo / |T”/ — Lf\dv'^ = 0, where has coefficients 


d{au * p„) -n_ d{bu * pn) 

d{v * Pn) ’ ’ d{v * Pn) ' 
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Indeed, recalling that ||/||(^i,2 < 1, we estimate 


/ — Lf dv"^ < J \a^{x) — a{x)\ dv'^ + J — 6| du'^ 

= J \iaiy*Pn){x)-a{x){u*Pn){x)\dx + I \{bu*Prr){x)-b{x){u*Pn){x)\dx 


< 2 


jUj{y - x)pn{y - x)dx 


iy{dy) = 2 / ijj{z)pn{z)dz 0. 


Thanks to this fact, we write 

limsup f It”/— T/| dz/"'= limsup 


£"/-£'7 




= limsup / |a"-a"| + |6"-6''n^ 

n—>-oo J ^ ' 

< J (|a — d| + |6 — 6|) dz^ 

where in the last step we apply (37). To conclude, it is sufficient to optimize upon a, 6, by 
density of continuous and compactly supported functions in T^(zz). 


A.4 Proof of Theorem 2.5 


We argue by iterating the three-steps scheme, the base case being that of diffusion operators 
with smooth and uniformly bounded coefficients. First, we extend the validity to uniformly 
bounded coefficients (without any regularity assumption), then to locally bounded coeffi¬ 
cients, and finally integrable coefficients. Although everything could be obtain in a single 
iteration, we think the approach highlights the different roles played by different approxima¬ 
tion procedures. Indeed, our crucial improvement with respect to [FigOS, Theorem 2.6] is to 
move from uniformly bounded to integrable coefficients, which is rather delicate: by compar¬ 
ison, in the deterministic case, one is able to deal directly with locally smooth coefficients 
(see e.g. [AGS08, Proposition 8.1.8]), essentially because paths either go to infinity, i.e., the 
solution explodes in a finite time, or stay in a compact set. Roughly speaking, the source of 
difficulties in the stochastic case is that we have to deal with “averages” of such behaviours, 
and moreover the solution to a genuinely stochastic martingale problem is expected to instan¬ 
taneously “diffuse” over compact sets (of course, with small probability as these sets become 
larger). 

Case of smooth and bounded coefficients. Let a, b be Borel maps as in (2), with 


l“i|lc2(Rrf) + \\bt\\c^(u<i) 


dt < oo. 


(44) 


Then, the superposition principle holds for every solution u = ^ .^(K'^) of the 

FPE (3). This follows from two well-known facts: existence of Ito’s stochastic differential 
equations and uniqueness for narrowly continuous solutions of FPE’s. 

The existence result is standard, with the possible exception of the integrable bounds 
with respect to the variable t G [0, T] (usually, one requires uniform bounds), but in fact such 
condition is sufficient for the various applications of Gronwall inequality. For the sole purpose 
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of establishing a case base for the superposition principle, the usual stronger assumptions on 
the coefficients, e.g. a, 6 G C'^((0,T) xR'^) would even be sufficient, at the price of introducing 
an extra mollihcation step with respect to the variable t G [0,T]. 

Theorem A.5. Let a, b be Borel maps as in (2), satisfying (44). Then, for every v G 
there exists a solution rj of the MP associated to L{a, b), with r]o = D. 

Proof. The assumption a G Ll{C^{W^)) entails that the symmetric non-negative square-root 
of a, i.e. the (essentially unique) map 

a : [0,T] X R*^ —)■ Sym^(R'^) such that = at, JSf^-a.e. t G (0,T), 

is bounded and Lipschitz with respect to x G R'^, with Lipschitz constant integrable w.r.t. 
t G (0,T), see e.g. [SV06, Lemma 3.2.3]. Then, it is sufficient to solve by Picard iteration the 
ltd stochastic differential equation 

dXt = bt{Xt)dt + at{Xt)dWt, Xo = X, 

where A is a r.v. independent of the d-dimensional Wiener process W. By Ito formula, the 
law of X, i.e. X^¥, is a solution of the martingale problem associated to T(a, 6). □ 

Of course, the MP is also well-posed, but we need a stronger uniqueness result, for narrowly 
continuous solutions of FPE’s, which is e.g. a consequence of results on backward Kolmogorov 
equations. We refer e.g. to the expository notes by [Kry99] for more details; notice however 
that, also in this case, the standard literature studies equations of the form 

dtf = -Ltf + g, in(0,T)xR^ h = f, (45) 

assuming a, b smooth and g G C'^((0,T) x R*^). A solution to the equation (45) is dehned as 
a function / G C'^’^((0,T) x R'^) such that 

dtf{s,x) = -Lsf{s,x)+g{s,x), for (s, x) G (0, T) x R'^, with lim/(s, x) =/(x). 

s\T 

To our purposes, we need existence of a solution, together the following regularity results for 
the solution / (which entails uniqueness): 

sup \\ft\\ci<{\\I\\cl+T\\g\\c2 )c( [ 

where z C{z) denotes some function depending on the dimension d only (the proof gives 
that C has an exponential behaviour). The proof follows by direct differentiation of the 
equation, see [SV06, Theorem 3.2.4] for a detailed derivation. Moreover, as a consequence of 
the maximum principle, if / > 0, and g > 0, then the solution / is non-negative as well. 

We are in a position to prove the following result, akin to [AGS08, Proposition 8.1.7]. 
Again, we provide a slightly stronger statement than what is needed for the superposition 
principle (e.g., we deduce uniqueness for possibly signed solutions of the FPE). 

Theorem A.6 . Let a, b be Borel maps as in (2) with 

+ \\bt\\c'2{B) every bounded open B C R'^, 

and V = (r't)tg[o,T] ^ .^(R'^) be a narrowly continuous solution of the FPE associated (3). 
7/t'o < 0, then vt < 0, for every t G [0,T]. Thus, for v G ./#(R'^) there exists at most one 
narrowly continuous solution v with uq = u. 



Iffilic? + ll^illcf 


dt + T\ 




(46) 
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Proof. Let g G C^{{0,T) x M*^), with g > 0: it is sufficient to show that J g dv < 0. Fix 
R> I large enough so that the support of g is contained in (0,T) x -B_r( 0) and let XR be a 
cut-off function, as below Remark 2.3. Notice that letting aji = axR and = bxR in place 
of a, 6, condition (44) holds and L^f = Lf on (0,T) x Rr( 0), for every / G Cl{{^,T) x R'^). 

For e > 0, let a|j, be a double mollification with respect to the space and time variables, 
and define = L(a^,6|j), which is a diffusion operator with smooth and bounded coeffi¬ 
cients, satisfying (44) uniformly in e > 0. Let be a solution to the backward Kolmogorov 
equation 

+ 5 ) It = 0 ) 

and choose f^XR in the weak formulation (5), which is admissible because G x 

Since /^ < 0 and z/q < 0; we have 


0 > - y fxR dvQ = j [xR dtf + £^ifxR)] dv 


= J [-Xr^rI + '‘^ifXR)]di' 

= J {xii [g + ^Rf - + r^XR + a(vr, Vxr)} dR 

> gdv- sup \\f!\\c 2 (Kd) / [XR\aR-a\ + \bR-b\ + \LxR\ + \a\\VxR\]d\R\. 

J te[o,r] J 

As e I 0, since or = a and bR = b on (0,T) x B{0,R), the second integral converges 
to /[|Lx_r| + |a| |Vx/j|] d |z^|, and sup^gjQ j’] ||/f||(^2 is uniformly bounded in e > 0, by (46). 
Finally, we let i? —)■ oo and conclude, since |Vxi?| -|- |Vxi{| —)■ 0, pointwise and uniformly 
bounded. □ 

The superposition principle follows immediately from these facts: since any weak solution 
R = (r't)zG(o,r) admits a narrowly continuous representative D, we let rj he a solution of the 
MP associated to L{a, b), with F = Fq (Theorem A.5) and notice that the curve g = {gt)t&[o,T] 
is a narrowly continuous solution of the FPE associated to L, with go = Dq. By Theorem A.6, 
we conclude that gt = I>t, for t G [0,T]. 

Case of bounded coefficients. We extend the validity of the superposition principle for 
diffusions with uniformly bounded coefficients: this already provides an extension of [FigOS, 
Theorem 2.6], as uniform bounds are imposed only with respect to x G R*^. Precisely, we 
assume that the coefficients a, b satisfy 


sup [|at(x)| -I- \bt{x)\] dt < oo. 

ajGK'* 


(47) 


Step 1 (approximation). We argue by convolution with a kernel p = aexp(—y^l -|- jxj^). For 
e G (0,1), let p^{x) = e"'p(x/e) and notice that IVV^^I < Ce~‘^p^, for i G {1,2}, where C is 
some absolute constant. Then, = r * p^ solves a FPE with respect to a diffusion operator 
with coefficients a^, P satisfying (the correspondent of) (44), as a consequence of the last 
statement in Lemma A.l. Existence of superposition solutions G ^{C{[0,T]-,R<^)) for the 
associated martingale problems follows from the smooth case settled above. 

Step 2 (tightness). We notice first that, being (i^o)e>o ^ narrowly convergent sequence of 
probability measures (thus, it is also tight), there exists some increasing function 0 : R —)■ R 
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with lim 2 _^oo = oo such that sup^^Q f 6{\x\)di/Q < 1. For R> 1, let then XR ■ [0,1] 

be the usual cut-off function and, for i G {1,... ,d}, let := XiXR £ ^■ For any (but 

fixed) p G (l,oo), let 0i(x) = © 2 ( 3 ^) = \xf, and apply Corollary A.4 to the solution rj^ 
(e > 0), with / := o a and these precise choices of 6, ©i, © 2 - We obtain some coercive 
functional 'h : (^([O, r];]R) ^ [0, 00 ] (depending upon 6 and p only) such that 

j '^{xRO-f)drj%x) < j 0{\x^R\)dr]^Q + + |af (Vrc*^, (48) 

Since ||Vx^||^ is uniformly bounded and is infinitesimal as i? —)• 00 , we may 

let R —7- 00 , and by lower-semicontinuity of 4', Fatou’s lemma and Lebesgue dominated 
convergence theorem, we obtain a similar bound with the functions re* in place of x^j^: 

I 'S(~i‘)dn‘(^) < I+ j(‘ (|{!>');|’’ + Ko'jJ’t) dvfdt, 

where we also make explicit the fact that Tj^ = vf, for t G [0, T]. 

Inequality (37) and the assumptions on 6 entails uniform bounds for e > 0 of the form 


J ^(7*)dr7=(7) <1 + J + l«rr ) dutdt. 

Tightness follows, for 7 1 —)■ Yli=i '^(7*) is coercive in C'([0,T]; 

Step 3 (limit). This step is fully covered in Section A.3. 

Case of locally bounded coefficients. Next, we assume that 

fT 

/ sup [|at(rc)| -|- | 6 ((x)|] dt < 00 , for every bounded borel B C 

Jo xeB 


(49) 


and we prove the validity of the superposition principle for every weak solution 1/ = {i^t)te(o,T) ^ 
.^(M*^) of the FPE (3) (recall that we also assume (4)). 

Step 1 (approximation). We approximate via push-forward by smooth maps. For M > 1, let 
Xm be the usual cut-off function and let ttm ; i-)- be the map 

= xxm{x), so that ttm{x) = x\m{x) G C'c(M‘^). 


By (49), it holds |T(7r)y^)| < ||vr]^^||c-2 sup|j.|<2M [|«(a;)l + \b{x)\] for x G f G {1,... d}, and 
similarly a(V 7 ^)^^, Vtt;^) < sup| 3 ,|< 2 M I«(a;)l 7 for x G M'', i,j G {l,...d}. 

Since conditional expectations reduce norms, we deduce that := {v) solves a FPE 

associated to a diffusion on M'^, whose coefficients a^, satisfy (47): thus the previous 
argument gives superposition solutions . 


Step 2 (tightness). The argument is very similar to the previous case, with the only caveat 
that ©1 and ©2 must be chosen more carefully. Indeed, we rely on the de la Vallee Poussin 
criterion, which improves the integral bound (4) to one of the form 



[©i(| 6 |) -F © 2 (|a|)] dutdt < 00 , 
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for some suitable ©i, ©2 that fulfil the assumptions of Theorem A.2 (the moderate growth 
assumption on ©2 can be always obtained, possibly passing to a slightly “worse” function). 
With such a choice of © 1 , ©2 (and building 9 as in the previous case, for {i'^)m>o is tight), 
we obtain for some coercive functional T such that inequalities akin to (48) and the following 
ones, with ©1(2:), ©2 (-2) in place of \z\^. 

Step 3 (limit). This step is described in Section A.3. 

General case. The final step consists of removing the assumption (49). 

Step 1 (approximation). We perform once again an approximation via convolution, e.g. as in 
the case of uniformly bounded coefficients. In this case, however, we only use the fact that 
(i/^)^ are solutions to FPE’s associated to diffusion operators whose coefficients are locally 
bounded (and the bound (4) is preserved). 

Step 2 (tightness). We argue exactly as in the previous case, i.e. using de la Vallee Poussin 
criterion to provide suitable © 1 , © 2 . 

Step 3 (limit). Again, this step is described in Section A.3. 

As already remarked at the beginning of this section, one could combine all the arguments 
above and prove Theorem 2.5, starting from the “base case” with a single combination of 
mollifications and push-forwards approximations. On a technical level, the main difficulty is 
to obtain the result for locally bounded coefficients, and this is done after we establish the 
result for uniformly bounded coefficients, regardless of their regularity, essentially because the 
push-forward approximation may not preserve it. 
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